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I. INTRODUCTION 


This research is concerned with a new realization of digital filters. The method 
used involves delta modulation. 

We define and discuss digital filters in Chapter Two. 

In Chapter Three we briefly explain delta modulation which converts an analog 
signal to a digital (binary) signal. We then introduce a new method of convolution 
which is called delta convolution in this report. The two signals to be convolved 
are first converted to binary form using delta modulation. Time domain analysis 
is used to derive equations of interest and to define the required hardware. An- 
alytical examples are presented to show that the results obtained using the delta 
convolution method and conventional methods are similar. We conclude that direct 
implementation of the delta convolution method results in a complicated system 
(large number of gates and interconnections). 

A ‘system which involves two feedback circuits to account for the previous 
values involved in delta modulation and convolution is developed in Chapter Four. 
We call this the reduced delta convolution (RDC) system. The result is a simple 
hardware realization. We define (block diagram) a system which functions as a 
convolver (linear filter) and a similar one which is a correlator. 

In Chapter Five we introduce a possible circuit to implement the RDC system. 
A comparison between the needed hardware and complexity of the RDC system 
and a conventional digital matched filter (DMF) is included. The result is that the 
RDC system has a simple layout for IC realization and requires less than 10 percent 


of the transistors needed for a typical DMF having comparable performance. 


18 


Simulation results when the RDC system is used as a convolver (lowpass filter) 
and a correlator (matched filter) are presented in Chapter Six. The simulation 
results include the transfer function of a RDC lowpass filter, the response of the 
RDC lowpass filter to two signals, one in the pass band and the other outside 
the pass band, and the output of a RDC lowpass filter with a square wave input. 
Included also are responses of matched filters to a raised cosine pulse and to a chirp 
signal. 7 

Chapter Seven documents the noise characteristics of the RDC system. Op- 
timum values of important RDC system parameters are derived. Analyses of self 
noise and input noise are presented. A conclusion made is that the RDC system 
has noise performance which is better than predicted by ordinary filter theory. 
Simulation results for the chirp matched filter and lowpass filter when used in the 
presence of additive input noise are shown. Again, the noise performance is better 
than predicted by theory. A reason for this is detailed. This chapter also includes a 
modification to the RDC method called variable step size RDC. Simulation results 
for the variable step size RDC lowpass filter and notch filter are presented. 

The results of this research indicate that the VLSI compatible RDC system 
may be an attractive means of realizing digital filters. The system requires no 
synchronization, has no multipliers and is easily programmable. Other features 
include real time operation (no off-line processing) and no limitations on the time 
duration of the signal to be represented in matched filter applications. 

When compared with other digital filter realizations, the RDC system appears 
to require considerably fewer components and simpler connections. Finally, because 
of the properties of delta modulation, the RDC system achieves noise performance 


which surpasses that predicted by ordinary filter theory. 
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II. BACKGROUND 


This research is concerned with the realization of linear systems (filters) using 
digital technology. Common realizations of digital filters require a large number 
ot same Busou and interconnections. This report describes a new type of we 
filter Saal documents its advantages and performance. We begin this descriptael 


by introducing linear systems. 


A. THEORY OF LINEAR FILTERS 
A continuous-time or analog signal may be described by a function of time, 


say z(t). The response of a linear time-invariant system to z(t) is given by 


y(t) = [ h(r) a(t —7r)-dr= [ x(\)- h(t — A)- dd (2.1) 


~od —oo 
where h(t) is the response of the system at time ¢ to a unit impulse applied at time 
t= 0. See Fig. 2.1. The function A(t) is, then, used to represent the linear system. 
Linear ore are used as convolvers to realize filters (lowpass, bandpass, notch) 
and as correlators (matched filters (MF)). 

The right side of Equation 2.1 is commonly called the convolution of two 


functions z(t) and A(t) and denoted as z(t) * A(t). From transform theory 


Ge ae | (2. 


IY) 
to 
SY 


where Y(f), X(f), and H(f) are respectively the Fourier transforms of y(t), x(t), 
and h(t). 
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A discrete-time signal is a sequence of numbers {z(z - T,)}, where the integer 
2 May vary over a finite or an infinite range, and where T, is the sampling interval. 
A linear time-invariant discrete-time system can be described by the input- 


output relationship 


ili SE ea aie ys Coe ni ((n—m)- T,) 
=T,- S> 2((n—k)-T,)-A(k-T,) (2.3) 
k=—0o 


where z(n-T,) and y(n-T,) are the input and output signals, and h(n- T,) is the 
impulse response of the system. That is, h(n-T,) is the response of the system at 
n-T, due to a unit sample (unit impulse sample) applied at t = 0. See Fig. 2.2. 

The right side of Equation 2.3 is the convolution sum of the two sequences 
{z(n-T,)} and {h(n-T,)}. When the sequence {h(n-T,)} has only a finite number 
of non-zero terms, we say the system has a finite impulse response (FIR). Otherwise, 
the system is said to possess an infinite impulse response (IIR). If h(n-T,) = 0 for 
n <0, we say that the eyeteal is causal or physically realizable. 

Using z-transforms 


ie ia) Xz) (2.4) 


where Y(z), H(z), and X(z) are respectively the z- transforms of y(n-T,), h(n-T,) 
pager (7 + 7, ). | 

A large class of linear time-invariant discrete-time systems can also be de- 
scribed by the linear constant coefficient difference equation 


N 


iF 
y(n-Ts) = So ax-y((n—k)-Ts) + D0 bea ((n—k)-T,) (2.5) 


hE k=0 
where it is possible to convert Equation (2.5) to an equation of the form of Equation 


(2.3) (Kirk and Strum [1]}). 





Fig. 2.1. Representation of linear time-invariant continuous-time system. 





Fig. 2.2. Representation of linear time-invariant discrete-time system. 


Dh 


In Equation (2.5), if a, = 0 for all k, then 
L 
Tome) = > bp aw((n =k))-T.) (2.6) 
k=0 
So, the output of the system depends only upon the current and previous values 
of the input. This system is called non-recursive. 
From these results, we can now consider ways of realizing discrete-time sys- 
tems. These realizations are called digital filters. This research is concerned with 


a new realization of a digital filter. 


B. VARIOUS REALIZATIONS OF LINEAR FILTERS 
The unit sample response for the non-recursive system described in Equation 


(2.6), is given by 


L 
h(n-T,) = Y_ dyd(n — k). (On) 


k=0 


This system is known as a finite-impulse response (FIR) system because we have 
a finite number of terms in the unit sample response (Bird [2]). 


From Equation (2.7), a FIR digital filter has a transfer function 


L 
ele ss bez” 
K==() 


which is often implemented as a tapped-delay line and is sometimes called a 
transversal filter. (See Fig. 2.3 and Rabiner (3]). | 

It is important to realize that, while a non-recursive filter must have a FIR, 
a FIR filter need not be implemented non-recursively. A recursive implementation 
of a FIR filter can be obtained by introducing poles in the transfer function and 
then canceling them with extra zeroes (Bird [2]). A new recursive implementation 
of a FIR filter, which is called the reduced delta convolution (RDC) method, is 


introduced in Chapter Four. 
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Fig. 2.3. Realization of transversal filter. 
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A transversal filter can have an arbitrary impulse response of finite time dura- 
tion and therefore can be used to implement any linear filter. (Any system whose 
output is linearly related to the input is a linear filter.) In this sense, a transversal 
filter can be thought of as a fundamental building block of linear systems. 


In Equation (2.5), if a, # 0 for one value at least, then we have an IIR filter. 


In classical signal processing most of the filters are IIR systems. These filters 
Poresaond to direct implementation of the filter’s rational polynomial enmetee 
function. 

In modern signal processing most filters are FIR systems. These filters corre- 
spond to a direct implementation of the filter impulse response approximated with 
a finite number of terms. These filters often arise as matched filters for radar and 
sonar signals (Kailath [4] and Rabiner [3]). 

FIR linear filters have several important properties which make them attrac- 
tive for digital signal-processing applications. Among these features are simple 
design, linear phase and the absence of any stability problems which may occur 
in IIR filters. On the other hand, for the same filter, non-recursive technique uses 
more taps and multipliers (Hamming [5], Young [6], McClellan [7], and Whalen 
[8}). 

Transversal filters are the popular choice of digital filter when using VLSI. We 


consider next various circuit realizations of transversal filters. 


C. TRANSVERSAL FILTER REALIZATION 


1. Block Diagram 


a. Transversal Filter as a Linear Filter 


To implement a linear filter using transversal filters, it is sufficient to 


choose the tap weights appropriately. From Equations (2.6) and (2.7), if the desired 
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impulse response is known, then the transversal filter tap weights are simply the 
sample values of the filter impulse response. See Fig. 2.4. 
b. Transversal Filter as a Matched Filter/Correlator 

Pulse detection in noise differs from analog signal transmission in two 
major respects. First, we are usually concerned with determining the presence 
or absence of a pulse. Second, we often know the pulse shape in advance, but 
not its Bade or arrival time. The pulse-shape information wales it posaiule to 
design optimum receiving filters for detecting pulses buried in noise having a known 
spectral density function. Such optimum filters (correlators) are called matched 
filters. 

Let z(t),0 < ¢ < Tq, be a signal to which a filter is to be matched. 


The impulse response of the matched filter is defined as (Turin [9] and Turin [10]) 
b-2(Tyg—t); O0<t<Ty 
h(t) = 
0; _ elsewhere 
where the gain factor 6 is arbitrary and is henceforth taken as unity and Ty is the 
pulse time duration. The name matched filter comes from the fact that h(t) has 
the same shape as the pulse z(t) with time reversed and shifted by Tz seconds. 
Fig. 2.5 is a diagram showing the application of a matched filter. 
The ratio of peak signal power S, to average noise power N, at the 


output of the MF is given by 


Noto! ae 


where S$; and N; are average input signal power and noise power and B is signal 


bandwidth. Then 
S»/No 


SNRI = 5,]N; = 





Breads (2.8) 


where SNRI is the signal-to-noise ratio improvement factor (Horrigan [11]). 
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Fig. 2.4. Diagram of a transversal filter as a linear filter. 
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Fig. 2.5. Diagram indicating the application of a matched filter correlator. 
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This important result shows that the filter SNRI is determined only by 


the time duration of the pulse and its bandwidth and not by the pulse ‘shape’. 


Fig. 2.6 shows a transversal filter implementation of a MF, where we can 


see that the values of the weights are exactly the same as that given for linear 


filters (Fig. 2.5) except that the positions of the weights are reversed. 


2. Circuit Realization 


In Figures 2.4 and 2.6 the delay needed for the analog signal can be im- 


plemented using the following conventional techniques: 


Charge-Transfer devices (CTD’s), which include both charge-coupled devices 
(CCD’s) [12 — 14] and bucket-brigade devices (BBD’s) {15, 16], can be used to 
delay the analog signals (Butler [17]). When CTD’s are used to delay analog 
signals, the signal to be delayed is first sampled at a rate greater than twice 
the largest significant frequency component of the signal. The analog samples 
are then clocked through the CTD shift register. CTD’s can achieve hundreds 
of milliseconds of delay (Buss [18]). 
Acoustic delay lines are used as alternatives to CTD’s for analog time delay, 
but for small values of time delays (< 20usec (Buss [18])). 
Digital shift registers (SR) preceded by analog/digital conversion (ADC) and 
followed by digital /analog conversion (DAC) can also be used to delay signals 
as shown in Fig. 2.7. | 

a. Transversal Filters Using CTD’s 


In order to make a CTD transversal filter, it 1s necessary to nonde- 


structively sample the delay line and to perform the weighted summation. The 


sampling is achieved in different ways depending on whether CCD’s or BBD’s are 


used (Buss [18]}). 


2g 





Fig. 2.6. Diagram of a transversal filter as a matched filter. 
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Fig. 2.7. Diagram showing realization of time delay of analog signals 
using digital circuits. 
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Generally, there are two approaches to implementing programmable 
transversal filters using CTD’s. 

The first is the analog/analog system (Bosshart [19], and Denyer [20]). 
Fig. 2.8 shows a block diagram of an analog/analog CTD transversal filter where 
analog stores for both input signal and impulse response (reference store) are used 
in conjunction with analog multipliers. 
‘The main limitation is that since the weighting coefficients are stored 
in the form of charge, they decay due to thermal leakage and must be refreshed 
every 10 to 100 ms. The tap weight resolution is limited by the fixed pattern 
noise associated with the MOS analog multipliers at each tap. In addition, the tap 
weight data that are stored in off-chip digital memory must be converted to an 
analog signal before it can be transferred to the on-chip analog reference store. 

The other technique is the digital/analog approach (Tiemann [21] and 
Tanak [22]) where the weighting coefficients are represented in digital form. See 
Fig. 2.9. This requires a multiplying digital-to-analog. converter (MDAC) at every 
filter point. The MDAC uses MOS transistors that route the output of floating 
gate ioe to either positive or negative summing buses. As a result of this, a fixed 
pattern noise results due to variations in the tapping transistor characteristics along 
the length of the array (Gandolfo [23]). 

We can summarize the main limitations of transversal filters realized 
using CTD’s as follows: 

e The time duration Ty of signals that can be processed using CTD matched 
filters is ultimately limited by the storage time of the devices (i.e., the time it 
takes a stored charge to be lost due to leakage). Typically the storage time is 
of the order of 1 second (Buss [18]). 
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Fig. 2.8. Diagram of the analog/analog method of implementing CTD 
transversal filter. (a) System. (b) System impulse response. 
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Fig. 2.9. Diagram of the digital/analog method of implementing a CTD 
transversal filter. 
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e The filter length (i.e., the number of delay stages L) is ultimately limited by 


charge transfer inefficiency. Calculations indicate that 
Oa ay (2.9) 


where L is the number of delay stages and ¢ is the charge transfer inefficiency 


per point where 


e = 107° for early versions of CCD/MOS 
and 


e = 0.4 x (10)~° for recent versions 


(Creasey [24}). 

e The Nyquist sampling theorem requires that a signal having a bandwidth B 
be sampled at a frequency greater than 2B. Combining this requirement with 
Equation (2.9) gives the following limitation on the Ty-B product (which is a 
measure of pulse compression ratio or processing gain) of signals that can be 


processed using CTD filters. 
ne baa L/e. (2.10) 


The limitation imposed by Equation (2.10) can be compensated for by selecting 
the weighting coefficients to invert the dispersion due to the imperfect charge 
transfer. But the dependence of charge transfer efficiency on the signal ampli- 


tude makes it impossible to exactly invert this dispersion at all signal levels. 
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e The signal bandwidth is limited to less than half the maximum clock frequency 

of the filter. 
o For CCD’s this limitation is + 20 MHz. 
o For BBD’s this limitation is + a few Megahertz. 

e Another limitation on the CTD transversal filter is the accuracy with which 
the weighting coefficients can be determined. For example the transversal 
Alters using MOSFET analog multipliers suffer from: 7 

o Poor accuracy 

o Lack of long-term stability 

o Drifts that are found in most MOSFET multipliers 
The weighting coefficients & also be realized using a split-electrode structure. 
However, because the weighting coefficients are fixed by the split-electrode 
structure, such devices are only suitable for applications where a fixed filter 
response is required (Haken [25]). 

b. Transversal Filters Using Digital Delay 

It is possible to use clocked shift registers to provide delay. Then in 
Fig. 2.10 the sampled input 2(n-T,) is quantized into a number of bits, typically 
12 to 16 bits for audio and 6 to 8 bits for video signals. 

The digital signals (code words) recirculate in virtual delay lines com- 
posed of random access memory (RAM) and are multiplied by the weights which are 
stored in programmable read only memory (PROM). Every sampling time period 
a new code word enters the RAM where the oldest one is discarded. Using a single 
multiplier, the multiplication rate has to be L- f,, where JL is the filter length and 
fs; is the sampling frequency. The large number of multiplications required each 


second is the main disadvantage of this scheme (Kailath [4]). 
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Fig. 2.11 shows an equivalent scheme which uses L multipliers and 
where m,; and m2 are the code word lengths of the PCM encoder and micropro- 
cessor system respectively. In this scheme the multiplication rate is reduced to 
f,. The main disadvantage of this scheme is the increased cost compared with an 
equivalent CTD implementation (Terrell [26]). 

The next chapter introduces a new approach to the realization of transversal 
filters. The circuit that results requires no voltage multipliers which result in 


a considerably reduced number of components (transistors) when compared with 


CTD and digital delay versions. 
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Il. DELTA CONVOLUTION 


A. INTRODUCTION 
This chapter introduces a new method of convolution which is called DELTA 


CONVOLUTION in this report. The two signals to be convolved are delta modulated. 


— 


B. DELTA MODULATION 

Delta modulation (DM) represents an analog waveform with a binary sequence 
(bit stream). This analog-to-digital conversion (ADC) is quite different from pulse 
code modulation (PCM). PCM involves code words and requires somewhat com- 
plex encoders and decoders. Recovering the analog waveform from the code words 
requires synchronization. By contrast, DM involves no code words, requires no 
timing signals, and uses simple hardware (Roden [27]). 

It has been found that analog signals such as speech and video signals gen- 
ell) have a considerable amount of redundancy; that is, there is a significant 
correlation between successive samples when these signals are sampled at a rate 
higher than the Nyquist rate. The redundancy in these analog signals makes it’ 
possible to predict a sample value from preceding sample values and to trans- 
mit the difference between the actual sample value and the predicted sample 
value estimated from the past samples. This results in a technique called DIF- 
FERENCE ENCODING. DM is one of the simplest forms of difference encoding 
(Ha [28]). 
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The idealized DM codec (coder/decoder) is shown in Fig. 3.1. The band 
limited output of the LPF is compared by subtraction with the stepwise approx- 
imation z,(t). The difference z, is passed through a hard limiter whose output 


1S 


ee). 


The resulting binary output Z(t) is applied through an ideal integrator (or accu- 


— 


mulator) with a feedback gain factor G to produce z,(t) where 








C= Az when using an accumulator 
i = when using an integrator; T is the sampling interval (3.1) 
and where 
ery [O25 2GHD): 2D) e441 
a ( 2 r, ~ when using an accumulator (3.2) 
= (iti Ge at when using an integrator 
where 


2n+1 





z4(t) = 24 ( T,) ;n-T, <t<(n+1)-7, Yn 


The decoder for a DM waveform is simply a staircase generator (accumulator) 
or an integrator. If a “1” is received the staircase increments positively. If a “OQ” 
is received the staircase increments negatively. The accumulator (or integrator) is 
followed by an amplifier having an amplification factor G. The amplifier may be 
followed by a low pass filter (LPF) to smooth the staircase output into a continuous 
function. A set of waveforms associated with the DM are shown in Fig. 3.2 which 
also defines Az of Equation (3.1). 

DM systems are subject to two types of quantization distortion which are gen- 


erally referred to as granular quantization noise and slope overload noise. Since 
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Fig. 3.1. Block diagram of the delta modulation codec. 
(a) Encoder. (b) Decoder. 
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Fig. 3.2. Delta modulation waveforms. 
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z,(t) cannot change by more than Az units in T; seconds, a is the highest input 
signal rate-of-change that the DM codec can follow. We call = the DM slope- 
following capacity and denote it by zg. When || exceeds this quantity, slope 
overload occurs and gives rise to the kind of error shown in Fig. 3.3. The granu- 
lar noise arises because the DM signal is a discrete-amplitude representation of a 
continuous amplitude process (Goodman [29] and Greenstein [30]). Consequently, 
granular eine is always present. Only granular noise occurs when the inoue 1s 
changing very slowly (* constant). 

The key to effective use of delta modulation is the intelligent choice of two 
parameters: (1) the step size Az, and (2) the sampling interval T,. Obviously, 
granular noise is reduced by decreasing Az, but at the expense of a reduced slope- 
following capacity, and, hence, greater slope overload noise. In Chapter 7 we 
introduce a procedure to choose the optimum step size and sampling interval such 


that the total noise obtained from the granular noise and the slope overload noise 


is minimized. 


~= 


C. DELTA CONVOLUTION 

The subject of this research is the representation of any linear filter using 
DM. Time domain analysis is used to derive equations of interest and to define 
the required fae In this chapter we pursue the DM equivalent of a linear 
filter by direct operation on the convolution integral. This results in a cumbersome 
expression and complicated hardware. However, working from this first result, it 
is possible to apply two successive modifications which result in a quite simple 


equation and a simple hardware structure. These modifications are presented in 


Chapter 4. 


4s 


Slope overload noise 
bandlimited analog input x(t) 


Staircase version x . (t) 
in AD 4 
slope = 


| 
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A x ato oo 





Esa peas granular noise 


Fig. 3.3. Waveforms showing slope overload and granular noise. 
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Assume that h(t) and z(t) are two analog signals where 
C3 Fey OE t<0 


Let A(t) and 2(t) be the delta modulated versions of h(t) and z(t), where h,(t) 
and z,(t) are the staircase representations of h(t) and z(t) such that we have no 


slope overload (Fig. 3.2). Assume also that 


— 


h; = i bit of A(t) = A(t) for i-T, <t <(i+1)T, (3:3) 
and 
x, = i bit of #(t) = Z(t) fori-T, <t <(2+1)-T, (3.4) 
We let h; and x; € {—1,+1} for all 1. Therefore, 
j 
h,(@)= Ah: She ee a (3.5) 
7=0 
and 
j 
a(t) Shey Saag j-1,<ta- Bee, (3.6) 
7=0 


When the sampling rate and step sizes are chosen properly, then, as seen in Fig. 3.2, 
cA (ts) ESF a Ita) j) Tet < 9 +1) (3am) 


Therefore, from Equations (3.6) and (3.7) we will have 
F ; 
z(t) o,(t) = Age) 23, pa Gee (3.8) 
7=0 
where | 
} = integer part of 
el r par —|. 
4 ger p T. 


In the same manner, 


h(t) m ha(t) = h- 7 hs 7 T, <t Gee (3.9) 


7=0 
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We now have expressions for z(t) and h(t) which involve the bits of #(t) and A(t). 


Now, a linear system (filter) is characterized by its impulse response h(t). 


Given an input voltage z(t) and impulse response A(t), the output voltage y(t) can 


be written as the convolution of x(t) and A(t). That is, 


Hey 1G) ed 0) / aan ener 


= OO 


— 


A conclusion is that in the time domain, linear filters are completely and uniquely 


identified by their impulse response A(t). 


The preceding development suggests that r(t) © z,(t) and A(t) = A,(t) and 


SO we can write 


YEAS hata, (ie [ ieee ab — 2) - dr 


—oo 


Let k = integer part of Ear Then from Equation (3.9) 


, | 
MG) An hn es kota <(k +1): 7, 


p=0 
Similarly, from Equation (3.8) 


z(t—r)=Az-) 23; m:-T,<t—r<(m+1)T, 


where m = integer part of a 





. But 7 = integer part of EA and so 


t=)7-734+ €1; €, is positive and less than Tj. 
T 


But k= integer part of Ea and so 


T=k-T,+ €9; €» 1s positive and less than T,. 
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(3.10) 


(3.11) 


(3.12) 


Gas) 


(3.14) 


From Equations (3.13) and (3.14) 


t — ae _ 
integer part of | = integer part of leuanenic masa 





Ts 
because €; and €, are positive and less than T’,. So, 
a i , 

the integer part of | 7 | = (3.15) 
From Equations (3.12) and (3.15) 

j—k 

a(t —G) AC mee. (GQ —k)-T, <t-—-r<(g—-—k+1)-T, (Sass 
7=0 


From Equations (3.10), (3.11) and (3.16) 


fore) k j—k 
ut)» | an Day ar | =| - dT 


p=—0 2=0 
Ly k jak 

~ / an ~~ hl ES -S a;| -dr (Sal) 
0 p=0 1=0 





where L, is the time duration of non-zero values of h(t). A duration of LZ, seconds 
corresponds tO 72 = = bits of A(t). 
A discrete version of y(t) is formed by replacing the integral in Equation (3.17) 


with another summation to obtain 


= to 7 <li (3.18) 


We write y(t) © y(j -T,) = Cy: { h(n) Aaj ie r)} - f= Oso. 


= 0 7 <0 
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where 


Cre Ar Nh hs 


and 


k=0 
where h, and zr, € {—1,+1}. 

We call h(r)Aa( j:T, —7) the delta convolution (normalized) at time oe Can at 
this last form, y(t) can be represented as the sum of appropriate products of +1. 
This result provided the initial motivation for this research because the product 
(+1) - (41) can be formed with XNOR gates. This is a considerable hardware 
advantage compared with usual digital multipliers (and usual digital filters) which 
have to form the product of PCM code words. 

As an example of delta convolution, we convolve h(t) with z(t) where A(t) and 
z(t) are given in Fig. 3.4. Applying the conventional methods of convolution the 
result y__ (t) is given in Fig. 3.5 when we let T, = 1. 

Assume that Az = Ah = 1. Then, the delta modulation versions of A(t) and 


z(t) are given in Fig. 3.6. Applying the delta.convolution Equation (3.18), 


y(0) = 0. 
The delta convolution when 7 = 1 (Fig. 3.7) is given by Equation (3.18) as 


TE Cae { a(n) A(T. 7 r)} 
— C7} : {ho : Lo } 
ney {1 1} 


= C, = 1 with our assumptions. 


49 





(b) 


Fig. 3.4. Waveforms. (a) Impulse response. (b) Input signal. 
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Fig. 3.6. Delta modulation versions. (a) Impulse response. (b) Input signal. 
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Fig. 3.7. The needed products for delta convolution at time T,. 
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The delta convolution when 7 = 2 (Fig. 3.8) is given by Equation (3.18) as, 
y(2T,) = Cr - {A(r)Aa(2T, — 7)} 
= C,- {ho(229 +21) +h 20} 
=C,-{1(2+1)+1(1)} 
ey 4) ee 
The delta convolution when j =3 (Fig. 3.9) is given by Equation (3.18) as 7 
y(3Ts) = Cr - {h(r)Aa(3T, - 7) 

= Cy: {ho(3z09 + 221 + 22) + hy (229 +21) + hero} 

= C,-{1(18+2+1)+1(2+1)4+1(1)} 

=O 

In the same way we can find the result of delta convolution y(t) for all values 
of 7 -T, (Fig. 3.10). It is clear from Figures 3.5 and 3.10 that simulated results of 
analog convolution and delta convolution are essentially the same for this example. 
Note that the maximum value of the convolution is 344.0 for both simulations. 
Theoretically, the maximum value of the convolution is 341.33. 

It is possible to expand Equation (3.18) and then create a circuit which forms 
y(j-T,) for any A(t) and Z(t). This is done in Appendix A. The result is a circuit 
having an unacceptably large number of XNOR gates and a large memory. We 
present in Chapter 4 a much simplified form of Equation (3.18) which eliminates 


the considerable redundancy of mathematical operations inherent in that equation. 


The practical result is a simple hardware realization of a circuit to form y(j - 


ayy 


o4 






Fig. 3.8. The needed products for delta convolution at time 27,. 
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Fig. 3.9. The needed products for delta convolution at time 3T,. 
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IV. REDUCED DELTA CONVOLUTION 


Implementing Equation (3.18) directly results in a large number of gates and 
interconnections as noted in Appendix A. The following reasoning suggests a sim- 
pler form may be possible. The delta modulation representation of a function (x(t) 
or A(t)) at any time t involves all previous values (history) of the function. The 
feedback portion of the delta modulator accounts for this history. Also, the result 
at any time ¢t of convolving two functions depends on the history of both functions. 
The integral portion of the convolution equation accounts for this history. 

When using Equation (3.18) to determine the delta convolution result at time 
t = t2,t3, etc., we essentially involve (calculate) repeatedly the previous history 
of both z(t) and A(t) at t2,t3, etc. This implies considerable redundancy of cal- 
culation (hardware) which increases with time. By carrying the history (previous 
calculations) forward. to the next calculation, the required hardware should be re- 
duced. Feedback can be used for this purpose. Indeed, use of one feedback circuit 
does a. 2 the complexity as shown in Appendix B. In this chapter, we develop 
from Equation (3.18) a system which involves two feedback circuits to account for 
the previous values involved in delta modulation and convolution. We call this the 
reduced delta convolution (RDC) system. 

The RDC system is a particularly simple implementation of Equation (3.18) 
involving only n, gates and two registers the lengths of which are n, and n; + 1 
stages where n, is the number of bits used to represent the system impulse response 
h(t). The form of the RDC system can be developed by considering the terms of 
Equation (3.18) at times T,,2T, and 3T, as follows: 


08 


y(Ts) = Cr - [20 - ho] Sr [Ts] 
lee dee = Ch : [ho : (229 is 2} ) ie hy : rq | 
=) 2 see 727, 


a y (Ps) ar E28) 


where 
al Ca = Cc; f [ho “21+ hy : ro| 
1 (2s) =7|7,) + r (27,] 
Ea po Ve 
Ht (bas) = ley 
and 
y (3T,) = Ci - [ho - (820 + 221 + 22) + hy - (2x9 +21) + he - zo] 
= 3r[T,] + 2r [2T,] + r [3T7,] 
= y (2T,) + I[3Ts] 
where 


r(3T7,] = C, - [ho - 2g thy -21 +ho- 29] 
1 3S ae ol ES | 
= I ([2T,] + r [8T;] 

Note that the present value of the output is the previous value plus an in- 
crement J. Further, the present value of J is the previous value plus a residual 
r which depends on the present contents of the input register. So, the new data 
updates r which then updates I which then updates the output. These updates are 


accomplished with feedback which involves storage of a discrete voltage as shown in 


Fig. 4.1 where y(j-7,) * y(t), and clock rate = 7 Although Fig. 4.1 has 
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—_ 


yj =1)i9 tj 


Stores read in on (k-1)th tick and read out on (k)th tick 


Fig. 4.1. Reduced delta convolution system. 
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been developed using transient or initial values (while the input register is becoming 
occupied), it is shown in Appendix B that this same form prevails when the input 
register is filled with values of z(jT,). So, the output of Fig. 4.1 represents the 
convolution of z(t) and A(t) when the input register is fully occupied, and the 
output is a transient response when the register is not fully occupied (loading and 
unloading). This transient behavior is part of any filter realization that uses delay 


— 


elements (digital or recursive filters). 


A. CIRCUIT REALIZATION OF TRANSVERSAL FILTERS 
USING THE RDC METHOD 


This section first considers a circuit based on Equation (B.11) which functions 
as a convolver (linear filter) and then a similar circuit which is a correlator. We 
show that a convolver or a correlator can be constructed by using two shift registers, 
ny XNOR gates, one or two amplifiers, one or two summers and two feedback loops. 
We show that by simply reversing the contents of the impulse response register, a 
convolver, in essence, becomes a correlator. 

1. RDC Transversal Filter as a Convolver 

From the preceding discussions and Appendix B, Fig. 4.2 is a diagram 

of a system which can be used to implement a RDC transversal filter which is a 
convolver. 

In the realization of the filter of Fig. 4.2, we show the necessary delays 

as consisting of charge-coupled devices (CCD). Not shown in Fig. 4.2 is the clock 

of rate T ticks/second used to transport the bit stream through register R, and 


used to transport the multilevel voltages through the CCD. 
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Fig. 4.2. RDC transversal filter convolver. 
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In Fig. 4.2, 


conv(j-T,) = output = y(j-T.) ¥ y(t) 
Aconv(j - T;) = conv(j -T,) — conv({y — 1] - Ts) 


A2conv(j - T,) & Aconv(j - T,) — Aconv({j — 1] - T,). 


In Appendix B, it is shown that 


— 


- nz,—l ny—l 
aweonn|)-T,| = Cy Ss een aren eae .(n) ) oe > | 
1=0 


—() 
where < Ry(z) > represents the contents (+1) of cell 2 of the impulse response reg- 
ister, etc.. and where < R,(z) > changes with time. So, the value of A2conv{j - T,] 
depends on the contents of register R, at time ) - T,. 


There are several constants in Fig. 4.2 which are defined as follows: 


Ny —] 


k; = amplifier voltage gain = —C, - SS De 
=O 


k 


amplifier voltage gain = Cy 
A;, Ao, A3 and Ay are voltage summers. 
2. RDC Transversal Filter as a Correlator. 
When the RDC system is used as a correlator, then the content of the 
impulse response register is reversed as shown in Chapter 2. Following the same 
procedures given in Appendices A and B, the following result is obtained 


nyol pel 
A2corr(j : ao) = =O ) » < R(t) Pane a UA? SF 1) ee ive 0) ae » | 


==) = 0 
where 
POCO) ee) 2 Acorr(j - T,) — Acorr( [7 —1]-T,) 


Meow!) secon) - faye corr |; — 1|- 7.) 
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and where A2corr(j-T,), Acorr(j-T,) and corr(j-T,) are the incremental change in 
the change of the discrete correlation values, the change in the discrete correlation 
values, and the result of correlation at time 7 - T;. 


Fig. 4.3 is a block diagram of a correlator as given in Equation (4.2) 


where 
Us| 1 


kz = C1 - S- h; = —k, and ky = —C, = —ko. 


=) 


— 


As expected, there is similarity of structure and commonality of hardware 
in the realizations of the convolver and correlator. In fact, if the input register is 
bidirectional and if the amplifiers become inverting, then a convolver becomes a 
correlator and vice versa. 

We can draw some conclusions concerning hardware from Fig. 4.2. 


a. When h(t) begins and ends at 0 volts, then* 


and we will not need the amplifier of aie k, nor the summer Ao. 

b. The needed storage equals (2n; +1) binary stages or cells and 2 single 
stage CCDs. 

c. The needed number of binary multipliers (XANOR gates) is n1. 

d. The needed number of operational amplifiers is 4. 

e. Cy 1s ascale factor (could be apart from the IC realization of the RDC 
system ). 

f. The same IC can be used as a convolver or a correlator. 

* Given h(n, -T,) =0. We know h,(n, -T,) & h(n, - T,). But 


ny—l 
h(m+T,)=Ah S| h; when 2(0) = 0. Sop whem h(r7pee— 0, then 
——') 
Risa! 
s h; = 0, under the condition h(0) = 0,(h,(0) is always 0). 


10 
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Fig. 4.3. RDC transversal filter correlator. 


We can now identify some important features of the RDC system. 

a. Uses digital technology (compatible with VLSI). 

b. There are no code words and hence no synchronization requirements. 

c. There are no multipliers (XNOR gates accomplish an equivalent function). 

d. Real time operation (no off-line processing). 

e. Easily programmable to represent impulse response of interest. 

ie | No hmitation on the time duration Tq of the signal in matched filter 
applications. 

g. No theoretical limitation on the number of register stages (taps); filter 
sections can be easily cascaded using several IC’s to realize longer 
duration impulse responses. 

h. Reduced hardware requirements (complexity) compared with other 
FIR filter realizations. An example is provided in Chapter 5. 

In Chapter 7, we conclude from simulation results that when used as a matched 
filter or ordinary linear filter, the RDC system has noise performance which is 
better than predicted by ordinary filter theory. The reasons for this are explored 
in Chapter 7. 

This chapter (4) presents a system (block diagram) which realizes an approx- 
imation to the convolution integral. The system can be used as a linear filter or 
a correlator. The next chapter details the hardware (circuit diagram) required to 


realize the RDC system. 
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V. HARDWARE REALIZATION OF THE RDC SYSTEM 


A. INTRODUCTION 

In this chapter we introduce a possible circuit to implement the RDC system. 
We also compare the needed hardware and complexity of the RDC system and 
a conventional digital matched filter (DMF). We conclude the RDC system has a 
simple layout for IC realization and requires less than 10% of the transistors needed 


for a typical DMF having comparable performance. 


B. CIRCUIT DIAGRAM OF THE RDC SYSTEM 
Fig. 5.1 shows a possible circuit having as output A2conv(j - T,) which is 
defined in Equation (B.11) and shown in Fig. 4.2. Register Rp, contains the bit 
sequence representing the delta modulation equivalent of the impulse response of 
the filter being realized. Register R, contains (n, +1) bits of the delta modulation 
representation of the input voltage to be filtered. The XNOR gates produce the 
binary products of the contents of register R, and register R,. From Equation 
(B.11), we need to sum the outputs of the XNOR gates. But, it is difficult to sum 
voltages. Therefore, we use field-effect transistors (FET’s) as voltage-to-current 
converters and then sum currents. Each FET will be in one of two states: 
1) Cut off when its gate has a “—1” (or 0) voltage which means zero current 
passing through it. 
2) Turned on when its gate has a “1” voltage which means I amperes passing 


through it. 
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Fig. 5.1. A circuit which provides A2conv(j - T;). 
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All the currents from the FET’s enter node B to form Jjo¢q; leaving node B. 

Bowirom Fig. 5.1, 

VY. = 13 + Leotal, 
or, 

Vi=r3-I- (na) 
where n,.& the number of bits for which < R,(7) >=< R,(7) >37=0,1,2,...ny —1 
Bey = n,— number of bits for which < R,(2) >#< R,(2) >; 7 =0,1,2,...n1 -1 
Down, —(n —Na)=2n,g—71; Ne =0,1,2,...n4. 


_As a result of this, 





4 = 2 (Atm); A= —n1,—7N, + 2,...,n, —2,n, (5.1) 
and, 
OSES Ses oe 
Peeume that (Vj)mar = 2V.-. Then V.. = rout where V,, is an externally 


supplied constant (dc) voltage. Using this in Equation (5.1) gives 


—— yee id 
) 2 





A+ Vac (5.2) 


Fig. 5.2 shows V; as a function of X. 
| We can use Equation (5.2) to define V2 in Fig. 5.1 as follows 
Hie i = f3°° I 

gene 








Va = Vee —Vi = Vor | 


Fig. 5.3 shows —V) as a function of X. 
ny—l 
Let k} = C; - ( = hs) = constant for the impulse response of interest. From 
7=0 
Fig. 5.1 


V3 = —V2 —ky- < R,(n1) > 


r3 





; Bs — ky- < R,(n1) = (5.3) 
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Fig. 5.2. Vi as a function of A. 





Fig. 5.3. —V2 as a function of . 
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From Equation (B.11) 


ny—l 


mviconv7-1,) =C, - yy ane) <A ?) >| eres (ti) 
pl) 
But, 
ni—l 
ee Oe 4-0) >=) 
7=0 


Therefore, we can write a 
Peon) 7 - Bel = Ciped — Ki < He(n).> - (5.4) 


From Equations (5.3) and (5.4), we can make V3 = A2conv(j - T;) if we set Cy = 
mae 
ae 
Knowing the step sizes Az and Ah, the sampling frequency = and the im- 
pulse response register contents and length n,, we can determine the constants C} 
and k,. Knowing n,, C, and choosing a convenient value for J, we can determine 
rz; and V,,. We choose Vpp in Fig. 5.1 sufficiently larger than (V} )maz so that the 
current J is constant. 
The circuit of Fig. 5.1 has as its output A2conv(j -T,). To find y(j - T,), 
the value of convolution, we have to use two feedback circuits as mentioned in 
Chapter 4. Fig. 5.4 shows one circuit which provides y(j - T;) when its input is 


A2conv(j -T,). The values of rs are arbitrary. From Fig. 5.4, 


Ul je ele=econ 7 = as) 


= conv ((j — 1)-T,) + Aconv() - Ts) 


a 
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where 


Aconv(j -T,) = Aconv((j — 1)T,) + A2conv(y - Ts) 
and A2conv(j - T;) is given in Equation (5.4). 


C. HARDWARE COMPARISON OF THE RDC SYSTEM AND A 
TYPICAL DIGITAL MATCHED FILTER (DMF) 


We can estimate the hardware requirements of the RDC system using Figures 
5.1 and 5.4. There are four different circuits identified in those figures. There is 
also a delta modulator which could be part of the system or could be off the chip 
when using IC technology. 

Assume 6 transistors are needed to realize each stage of a shift register. Then 
the RDC circuit requires 6n; transistors for the impulse response register and 
6(n, + 1) for the input register for a total of 6(2n, + 1) transistors. 

Assume each XNOR gate requires 10 transistors. Then the array of gates 
requires 10n, transistors. 

There are n, field-effect transistors. 

Assume 20 transistors are required for each op-amp for a total of 120 transistors. 

The total hardware requirements, neglecting resistors and capacitors, is (23n;+ 
126) transistors plus 2 single stage analog storage elements (CCD’s) and possibly 
a delta modulator. The hardware (in fact, the circuitry) is the same for the re- 


alization of a convolver (lowpass filter, bandpass filter, etc.) as for a correlator 


(matched filter). 
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A conventional way to implement a DMF is presented by Turin [10] (see Ap- 
pendix C). We will compare the hardware required for this DMF with that previ- 
ously tabulated for the RDC system. 

Assume that the sampling frequencies for the conventional method given by 
Turin [10] and the RDC method are 2fmaz and 2-F'- fmaz respectively, where 
fmaz 18 the maximum frequency component of the band limited version of the 


- — 


input signal and F is the oversampling ratio. As a result of this 


where 7n, is the register length of the RDC system and L is that of the conventional 
DMF register. 

We can now list the needed hardware to implement a conventional DMF and 
that needed to implement an equivalent RDC filter. The results are in Table 
.1. 


If we assume each. adder needs only 6 transistors, then the needed transistors 


~ 


to implement the DMF filter, neglecting the D/A and A/D converters and the 


amplifiers required are 
6(N+M)-L+10M-N-L+6M-N(L-1) 
where the needed transistors for the RDC system are 
23n, +126 = 23(F'- L) + 126 
Define Rr as 
the needed transistors when we use RDC method 


Rr 


~ the needed transistors when we use DMF method 
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TABLE 5.1. COMPARISON BETWEEN THE DMF CONVENTIONAL 
METHOD AND THE RDC METHOD 


IN 


Conventional DMF 


N shift registers having 
N - ZL binary stages and one 


having M.-L binary stages 


M-N-L binary multipliers 
(XNOR gates) 


M.-N-(£—1) summers 


A/D converter 
D/A converter 


M -N sets of interconnections 


RDC 
Two shift registers having 
2n,; + 1 binary stages 
n, binary multipliers 
(XNOR gates) 
6 operational amplifiers, 
each having less than 20 
transistors. 
Delta modulator 


none required 


n, interconnections 


where the number of these sets 
equals to L 


For large values of L, then 


23F 
tr ® oy 4 M)+ 16M N (5:5) 
A common value for M and N is 8 
235e | 
th ~~ — 
en Rr 1120 (5.6) 


From Equation (5.6) we conclude that the circuit complexity is equal for F' = 50. 
When L = 1000,M = N = 8, and F = 4 then Rr = 0.0823. In general, we 


conclude that the number of transistors needed to implement the RDC method 


19 


is less than 10% of that needed to implement the conventional DMF filter. In 
addition the circuit layout of a conventional DMF method is more complex than 


that of the RDC system. 
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VI. PERFORMANCE OF THE RDC METHOD 


In this chapter we present simulation results which show the transfer function 
of a RDC lowpass filter (LPF), the response of a RDC LPF to two signals, one in 
the pass band and the other outside the pass band, and the output of a RDC LPF 
with a square wave input. Included are responses of filters matched to a paige 


cosine pulse and to a chirp signal. 


A. LOWPASS FILTER PERFORMANCE 


1. Transfer Function of the RDC System 
Fig. 6.1 is a block diagram of the RDC system used to find the transfer 
function when used as a LPF. We assume a LPF cutoff frequency f, of 1.0 Hz. 


The impulse response A(t) is, then, 


sin(27f,t) sin(2zt) 


Wt) = 2 fo Oe Ot 


where the corresponding causal impulse response h(t) together with its staircase 
representation h,,(t) are plotted in Fig. 6.2. 


We obtain the transfer function by applying to the RDC system the signal 
Ee ea, cin 7f; t). 


Following the transient time, which is dependent upon the register length 


n, and the sampling interval T,, the output y(t) is 
y(t) = Ao - sin(27f;t + ¢) where ¢ ~ 0. 


The ratio 40 as a function of frequency is taken as the magnitude of the RDC LPF 


transfer function. 
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output 






analog 


input / 


x(t) =A -sin( w-t ) 






delta modulator 








Fig. 6.1. Block diagram of the system 
used to find the RDC LPF transfer function. 
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The register length n, depends on T, and Ty which is shown in Fig. 6.2. In 
facta a. For a fixed T,, increasing n; permits a more inclusive representation 
of the impulse response h(t). In the following, we measure Ty in terms of the number 
of zeroes about the main lobe of h,(t) and included in Ty. 

Figures 6.3, 6.4, and 6.5 show the magnitude of the transfer function for 
the RDC LPF where the number of zeroes of the impulse response are 2, 6, and 8 
respectively. As expected, the LPF transfer function “improves” as more zeroes of 
h(t) are included. 

Fig. 6.6 shows the magnitude of the transfer function of the RDC LPF 
when Tq includes 6 zeroes of h,(t) and when rectangular and Hamming windows 
are used. We see that the ripple in the transfer function is reduced using the 
Hamming window. 

When the step size Ah becomes larger than the side lobes of the impulse 
response, then the delta modulator considers the side lobes are of equal amplitudes 
(see Fig. 6.2). To alleviate this problem we can choose a smaller step size and 
increase the sampling rate s maintain the DM slope-following capacity. This 
results i increased register length n,. An alternative is use of a variable step size 


RDC system which is treated in Chapter 7. 
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2. The Output of the RDC LPF When the Input is the Sum 
of Two Sinusoids, One in the Pass Band and the Other 
Outside the Pass Band 


Assume that the RDC LPF with f, = 1.0 Hz, has an input 


t 
z(t) = 0.5 sin (= + 1.0 sin(w,t) 


where w, = 27f; and f, = 1.2Hz. 


Figures 6.7 and 6.8 show the input z(t), its staircase representation z,(t), and the 
output y(t) of the filter. We see that the 0.4 Hz sinusoid is preserved in the output 
while the 1.2 Hz term is suppressed. The RDC system when configured as a LPF 
does discriminate (filter) on the basis of frequency. 
3. The Output of a RDC LPF with a Square Wave Input 

Assume f, = 1.0 Hz. Assume the input z(t) is a square wave of frequency 
fo = 0.25 Hz, where Fig. 6.9 shows the input z(t) and its staircase representation 
z(t). It is expected that 7 first and third harmonics of the input will appear in 
the output of the LPF and all other harmonics will be suppressed. Further, the 
amplitude of the third harmonic should be one third that of the first harmonic 
according to the Fourier series representation of oe square wave. The simulation 
result agrees with this expectation as shown in Fig. 6.10. 

Assume that the square wave input has a fundamental frequency of 0.5 
Hz. Fig. 6.11 shows that all harmonics are suppressed in the output of the RDC 
LPF, as desired. 
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B. MATCHED FILTER PERFORMANCE 
1. Filter Matched to a Raised Cosine Pulse 
Fig. 6.12 shows a block diagram of the system used to simulate the im- 
plementation of a filter matched to a raised cosine signal. The impulse response of 


the matched filter 
A (1+ cos (20 (2 - 3))) Stes 


A (1+cos (20 (32 +4))) = SES, 


d 


— 


h(t) i CE — 5 = 


where z(t) is the input signal and Ty is the time duration. Fig. 6.13 shows the 
raised cosine impulse response where A was chosen to be equal to 1.0 volt and 
eye) 1.0 second. 

Figures 6.14 and 6.15 show the results of correlation when we use the 
conventional methods [y_, (t)] and the RDC method [y(t)] with an oversampling 
ratio fF = 8.0. Note that these results have the same shape. The theoretical 


maximum of the correlation is 


_ = APT, = 1.5 


Simulation results give a maximum value of correlation equal to 1.49999 when using 


the conventional method and 1.54536 using the RDC method. 
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raised cosine signal ; RDC output 
delta modulator raised cosine 
MF y(t) 





x(t) 


Fig. 6.12. Block diagram of the RDC matched filter that was simulated. 
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2. Chirp Matched Filter 
The chirp or linearly frequency modulated signal has been popular in pulse 


compression applications in the past because it is possible to realize filters matched 
to the chirp signal using acoustic devices. The instantaneous frequency f;(t) as a 
function of time is given in Fig. 6.16 where f;(t) = fi + arly where Ty is the 


chirp pulse time duration. So, the chirp signal (voltage) is then equal 


rn . 


Ase (wit i eee) 


A block diagram of the system simulated is shown as Fig. 6.17. Figures 
6.18, 6.19, and 6.20 show the input signal x(t), the output y.,,(t) when we use 
a conventional method, and the output y(t) of the RDC matched filter. System 
parameters used are f/f; = 0.0 Hz, f. = 5.0 Hz, Ty = 4.0 seconds, A = 1.0, F = 4, 
and n; = 160. 

From Figures 6.19 and 6.20 we observe that the maximum of y,,,(¢) is 
1.88820 and of y(t) is 1.84669 where the theoretical maximum is 1.8879557. After 
the input expires the output is constant (non-zero). This dc offset is explained 
in Appendix D. From Equation (D.3) the dc offset is expected to be equal to 
+42 (total area under h(t)). 

Increasing the oversampling ratio Ff’ to 32 decreases the dc offset. See 
Figures 6.20 and 6.21. This is so because when F’ increases, Ax decreases to keep 
the value of the slope-following capacity constant. And when Az decreases, then 
from Equation (D.3), the dc offset decreases. 

Note that in ine case of the raised cosine matched filter, the result of 
correlation has no dc offset because the total area under h(t) equals 0 (see Fig. 


6.15). 


96 





Fig. 6.16. The instantaneous frequency of the chirp signal. 


chirp signal 






X(t) 


Fig. 6.17. A block diagram of the RDC chirp matched filter that was simulated. 
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VII. NOISE PERFORMANCE OF THE RDC SYSTEM 


A. INTRODUCTION 

A DM converts an analog signal to a binary signal. This implies quantization. 
As a result of this and because of noise which may appear on the input, the output 
is subject to distortion. We identify two types of distortion and call their cannes 
self noise and external noise. Self noise is caused by the delta modulation process. 
External noise is present at the input to the DM. 

In this chapter we consider the effects on the output of the RDC system of 
first self noise and then external noise. Simulation results for lowpass filters and 
matched filters are presented. The idea of filtering on the basis of the slope of a 
voltage is introduced. Also included is the variable step size RDC method used 
to overcome a self noise effect on the representation of the impulse response of a 
filter. 

Self noise is affected by the step size and sampling ee used in DM. Therefore, 


before considering self noise, we first establish an optimum step size. 


B. OPTIMUM STEP SIZE 

In this section we define and derive the optimum step size for use with the 
RDC method. 

Let z(t) be a zero-mean random input to the two level quantizer shown in Fig. 
7.1. (In Fig. 3.1, the hard limiter, amplifier of gain G, and the flip flop can be. 


replaced by a quantizer in the sense given in Fig. 7.1). 


OZ 


QUANTIZER 
x y = Q(x) 





(C) 


Fig. 7.1. Two level quantizer. (a) Block diagram. (b) Transfer characteristics. 
(c) Example of input pdf. 
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Let z(t) have a probability density function (pdf) p,(-) and variance o?. (If 
the input z(t) has a nonzero mean, we can subtract it from the input and add it 


back after quantization.). The quantizer output y is 
+Az if x € (0,00) 
y= Q(z) = (7.1) 
—Az if x € (—o0,0) 
The quantization error ge = © — y is a random variable with pdf p,(-) and 


variance a7, where 


a, = E {42 | qo - Pq(q)- dg 


But, 


Ge = x — Q(z) is a function of z. 


So, 


7 = Bae] = | *(e— Q(2))? -pa(2) de 


[e —(-Az))* - p(x) - dz + [¢ — Az)’ - pz(z)- dz 


— =—o2 an 0 
- pi) aa r+ p,(r)- dz 
+ (Az)? - / p2(z)- dz —2Az- / Espa aa 


os 0 
If p,(z) is symmetric about zero, then 


OO 


oe =o*—4Az- / z-p,(z)-dz+(Az)’. (Tea) 
0 : 


We define the optimum step size Ar,,; as that which minimizes Oe To find Azop:, 


form 
do? 


q _ — ° ° —_= 
(Az) ~ | © p2(2) dr-t 24g — 6 
0 
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and solve for 


Pape — 7 De oan. (73) 
0 


Using this value in Equation (7.2) gives 


Te = 02 — 2Az-Arzoy: + (Az)? when Az # Azop: (7.4) 
and : = 
min {7} =a (Axon) when Az = AZop: (7.5) 
If we define a quantizer performance factor ¢? as 
2 
o 
e= 5 (7.6) 
then, 
. 2 ») 
. 2) _ min{aj} (Azopt) 


Y z 


where the inverse of =; can be considered as a signal to noise ratio. 

Table 7.1 lists Sues for different probability density functions and indicates 
the corresponding minimum value of Exe When the input is a sine wave, then the 
pdf is an arcsin function and 4222 = 0.9 (Jayant [31]). In all simulations this 
value of optimum step size was calculated for the applied voltage and then used to 
obtain the results. ; 

The value of Az,,; is that which minimizes the variance of the difference be- 
tween the input and the output of the two-level quantizer. The two-level quantizer 
is part of a delta modulator. In fact the DM of Fig. 3.1 is formed by an integrator 


or accumulator in a feedback path from the output of a two-level quantizer to its 


input. We now have a new variable which is the gain g in the feedback loop. In 





the next section, we derive the optimum value of g. 
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TABLE 7.1. OPTIMUM STEP SIZE AND MINIMUM VALUE 
OF THE QUANTIZER PERFORMANCE FACTOR FOR 
DIFFERENT INPUT PDE’S. 


pdf a min{es} 
Uniform at 0.25 
Gaussian J2 0.363 
ite 
Laplace Fr 0.5 
Gamma wa 0.667 





1. Optimum Gain of DM Feedback Loo 





In Fig. 7.2, assume that the accumulator has an inherent 1 bit delay; then, 
eT, gee ee 
As A result of this 
| Oe Te pay) cg Aare de =2z(n-T,)-—g-z({[n—- 1] - Ts) 
(2, (7-F,))" (ee eg a(n) cre 


Therefore, 


2 eA 2 « 2 Font 
=G, 19 “0, —29 ioe (7.8) 


DY 


Oo 


where o% is the variance of the input signal to the quantizer (hard limiter) and p; 
is the value of the variance normalized autocorrelation function R,,(T) of x(t) for 


7 = i 


Pe 
Setting ss = 0 yields the optimum value of g and min{o4} the minimum error 


variance. 
f ] 
Oon 


Og 


= (2g — 2p1)- 0, =0 
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delta modulator 


fl, 







accumulator 


e 
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delta demodulator 


accumulator 





channel 





Fig. 7.2. Block diagram of the delta modulator and demodulator. 
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Therefore, 
Jopt = pi and min{or} = (1 — pj)oz (7.9) 


We see from Equation 7.9 that the error variance is less than the input signal 
2 
variance a? for all non-zero values of p;. We call = a prediction gain Gp. 


2 
os 


o- 


The maximum value of Gp = mazr{Gp} = = (1 — p?)7? C7} 


min{a% } ~ : 

Processors of speech and images are often based on long-term autocorre- 
lation functions. In general, for speech p; has a value close to 0.9 (Jayant [31]). In 
sub-optimal design, p; is taken as 1 and so g = 1. 

The preceding consideration of a two-level quantizer and feedback gain 
can now be applied to determine the optimum step size for DM analog-to-binary 
conversion. 

2. Optimum Step Size for Delta Modulation 
From Fig. 7.2 we can define for a DM signal the reconatmuction error 


r(n- T,.) as follows 

ss T(t eg — De re) 
where y(n -T,) is the output of the delta demodulator which is equal to z,(n- T;). 
We can prove that the reconstruction error has a minimum variance 


min {€j}-(1—- pi) 
(1 — pj-min {€2}) * 


where we note that the minimization in Equation (7.11) is with respect to Az and 


min {a2} = ou Kora = 


(7.11) 


that this result holds for g ¥ p; (Jayant [31)). 


From Equation (7.7) we can say that 


ALopt = (\/1 = min {3} “On 


anon (7.12) 
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where © is a constant which depends on the pdf of the input z(t). Also, from 
Equations 7.7 and 7.11 


a pi) 2 


= pf min (EB])™* pe 


i 


So, from Equations (7.12) and (7.13) 


a | en, ae | 
Azopt = © (1 — p?- min {e2}) ” (7.14) 


where min{é€7} < 1, and can be found using Equations (7.3) and (7.7) or Table 
7.1 for the pdf of interest. | 
If p} — 1 (which is the case for speech and video signals for sampling 


frequencies higher than the Nyquist frequency), Equation (7.14) can be rewritten 


a ey, ,| a 
AZ opt © (1 — min{€?}) bo 


Then using Equation (7.12), 


: Sept an) ey Ox: (15) 


For pdf’s of many common input signals, 


as follows 


«x F? 





mazr{Gp} = 
ia, =i 


(Jayant (31]) where F is the DM oversampling ratio defined in Section V.C. 
Let 
mamiGp) — C4 = : (7.16) 


1 — p? 


where C,, is a constant. From Equations (7.15) and (7.16) 





AT opt = Cia) 


For example, for a sinusoid of frequency f; and amplitude A which is sampled at 


frequency f, 
r(n-T,) = A-sin(w,-n-T;) 





Therefore, 
R,,(T,) = A*sin(w; -n-T,) -sin(w; - (n+ 1] - T,) 
vee 2 , 
= cos -T,) — cos(w1(2n + 1) - T;) 
2 
= — cos(w, - T,) 
= 005 (25ft) = © cos (=) 
= oe 7, J = 3 cos { = 
and 
A? 
fee 0 os 
So 
Ross) TT 
= Teor © eo 
But 
x re . rt 
cos rz = 5 +o 
Therefore, 
| liny as 
Pi = cos @ ~1- 5 e. where f >> 1 
5G 


As a result of this 


From Equation 7.16 


1 
Cas = es 


(7.20) 





From Equations (7.17) and (7.20) 


T 
Nine — pate (e213) 
But 
A 
z /2 : 
Therefore, : 


Having established optimum values of step size and gain for a DM, we now 


can use these results to determine values of self noise in the ADC process using 


DM. 


C. SELF NOISE 
Chapter Three identifies two types of self noise: granular (quantization) noise 
and slope overload noise. 
1.. Granular Noise Power 
If |x(t)| is less than the step size as shown in Fig. 7.3, the granular noise 


is given by 


é, — 20) rd) 


Generally, we have granular noise for all x(t) because of quantization. To 
precisely define e,, let us assume that the analog input to the DM modulator z(t) 
and the output of the DM receiver are given in Fig. 7.4. Assume also that the 
step size Az and the sampling interval T, are infinitesimal, but with their ratio the 
same as the actual system. The output of the DM receiver will then be the smooth 


function z,,(t) shown in Fig. 7.5. 


ak 


We define slope overload noise €,, as the difference between z(t) and z,,(t); 
the remaining distortion (z,,(t) — z,(t)), 1s granular noise. 

Granular noise €, is essentially uncorrelated with the signal (Peebles [32] 
and Protonotarios [33]) and is approximately uniformly distributed (Goodman [29] 


and Protonotarios [33]), over the interval (—Az, Ar). See Fig. 7.6. 


Agr 1 = 
Blea) = fee (sa) — 


—-Azr 


Therefore, 


and the granular noise variance 1s given by 
Ag 2 
1 (AG) 
ue = / (€,)° (sz) dé, = 3 : (72235) 
—Ar 


So, the granular noise power increases as the square of the step size. 





2. Slope Overload Noise 
In general the slope overload noise is correlated with the message (Peebles 
[32]). A number of analyses to find the slope overload noise power are introduced 
by Protonotarios [33], Zetterberg [34], Rice (with O’Neal) [35], and Abate [36]. All 
these results do not either individually or collectively pertain to all slope-following 
capacities and input spectra. Greenstein [30] found an expression for slope overload 
noise power that is accurate for all slope-following capacities and input spectra of 
possible interest for linear delta modulators having Gaussian random inputs. 

Let us define the slope overload factor 
34 (=) (7.24) 

dt/r.m.s 
where z! is the slope-following capacity. A highly accurate approximation to the 


slope overload noise power o%, for linear delta modulators having Gaussian inputs 


and for all spectra, that was given by Greenstein [30], is as follows 
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Fig. 7.3. Illustration of the case where |z(t)| < Az. 


113 





Fig. 7.4 The input signal and its staircase representation. 


x(t) 





Fig. 7.5. The input signal and the smooth staircase representation. 
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Fig. 7.6. Probability density function of granular noise. 
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o?, = 0% [1 + 2.7538 + 2.95287] exp(—0.34157) 


3 


-exp {(a; — 2.753)s + a2 [exp(az3s + aqs”)-—1]};0< 5s <4.0 (7.25) 


2 


~, with an accuracy of 


where the expression given by Equation (7.25) predicts o 
1 dB or better for s between 4.0 and 6.5. For s > 6.5, 07, is at least 119 dB below 
o*. The variables a1,a2,a3, and a, are determined by the spectrum of the input 


(Greenstein (30]). For example, for bandlimited white noise 
a, = —0.036, a2 = 0.37,a3 = —3.83, and ag = —5.9. 


(See Fig. 7.7.) From this section and section VII.C.1 we can draw an error variance 


curve for DM (See Fig. 7.8.), where 


oi =0%, +.0% (7.26) 
Note that as Az falls below the value Azop;, the distortion increases more rapidly 
than when Az takes on values greater than Azo,:. We say that the delta modulator 


“fits” the signal if Az = Azop:. 


3. Signal to Noise Ratio Improvement (SNRI) Due to Delta 
Modulation 


Assume that we have chosen the step size to be equal to Az,,; where we 


have as input of the RDC system the signal alone. Then, 


Z 
2 wo} = tort) (7.27) 


where o2 ,07, , and o?_ are the total quantization noise power, the slope overload 
z x Jz 


noise power and the granular noise power respectively. 
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If we have as input noise alone (which has power a2) and the step size is 


=O. ie (7.28) 
where a? ,o3,_, and ae are the total quantization noise power, the slope overload 


noise power, and the granular noise power respectively. 
Neglecting the granular noise, which is practically very small compared to 
the slope overload noise, the effective noise to the RDC system will be decreased 


by o%,,- That is, 


the effective input noise power = 07, — 05, = 07(1— A) (7.29) 
where 
Ts0n 


See Fig. 7.9. From Equations (7.25), (7.29), and (7.30) we can find the input 
effective noise power when we have noise as input. 

When we have as input to the RDC system the signal as well as noise, the 
root mean square value of the slope of the input will increase. From Equation (7.24) 
as ()_. increases (where Az and T, are not changed), s decreases. From Fig. 
7.7 as s decreases, the slope overload noise power, in this case o2,_ ,,> lcreases. 


Neglecting the granular noise, the output power of the RDC system (neglecting 


any SNRI due to matched filtering) is 


2 


2 2 
or 50 on — 7 30240 


Therefore, 


eS (7.31) 
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where 


o2 
D= = ai 


mr 


So, the signal to noise ratio improvement due to using delta modulation = 5 pb 


D. THE ERROR IN CONVOLUTION WHEN THE OUTPUT BITS 
OF THE DELTA MODULATOR ARE IN ERROR 


The’DM is an ADC. The bit stream output can be transmitted as a baseband 
signal (twisted pair, coaxial cable, or fiber) or as a bandpass signal (modulated 
carrier). Errors may occur in transmission. 

A known signal at the input to a DM creates a known bit stream. Noise added 
to the signal creates a different bit stream. The difference can be considered as bits 
in error. The relationship between the noise characteristics and the bits in error is 
not generally known. 

In this section, we consider the effect on the output of the RDC system of bits 
in error. Because there are two feedback circuits in the RDC system, bits in error 
on the input affect future values of the output. 

From first principles of delta convolution (Appendix A), the error in the con- 


volution due to one bit in error is given by 


ny—l 
Sy os oo OF \- hy: (ny - o| ie) 
1=0 
= +€, 


where €, is a constant for a given A(t) and T;. 
Assume that we have two bits in error. Then the bits in error can be as follows. 


Both can be +1; both can be —1; or they can be of opposite sign. There are 


Z 
( ;) combinations of two 1's 


i 


2 
( 4 combinations of a single 1 and a single — 1, 


and 


@ combinations of two — 1's. 


where (™) is the binomial coefficient of m things taken n at a time. 
n 


Therefore the error in convolution will be 


: (2) 


2€, with probability on 


eID le 


2 
0 with probability G = 


and 


(0) _ 


—2e, with probability 52 


Ht |e 


Generally, if we have m bits in error and m is even, the error in convolution 


is given by 
m+ €1 with probability Ca) 
(m—2)-e€, — with probability aa, 
(m—4)-«€, | with probability oe 


De with probability = 
0 with probability a 


(m — k)ey with probability — 


and —m-e, with probability 
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where k is even number. See Fig. 7.10. 


So, 


Econy = 0 


and 





Fn = 2. (m— 8-4) tat (7.33) 


| Sees () 2 


~ —~— 


where m and k are even numbers. 
In the same way, if we have m bits in error, where m is an odd number, the 


error in convolution is given by 


ee) 


(m — k)-€, with probability oe 


So, 


E€conv = 0 
and 


= FY (Om = 8) 0) (." ) | (7.34) 


where m is odd and k is even (see Fig. 7.11.). Note that Equations (7.33) and 
(7.34) are exactly the same except that m is even in Equation (7.33), and it is 
odd in Equation (7.34). Note also that the error in convolution, for a certain error 
pattern (where all bits in error have left the first n, stages of R,), is independent 


of time. 


When the bits in error propagate through the first n; stages of register fz, 
the error in pence lation is dependent upon the position of the bits in error and so 
it is dependent on time. We want to find the signal to noise ratio at those times 
when the bits in error have passed through register R,;. Assume that we have m 


bits in error. Therefore, from Appendix A, 
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Fig. 7.10. The probability of error in convolution, m even. 
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p ( poate, 





Fig. 7.11. The probability of error in convolution, m odd. 
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ni—l pas 
the output + C; - | yy hi(ny — } | Se ¢\-- nose at (= 77 ae 
1=0 


1=0 





where the noise has zero mean and variance given by Equations (7.33) and (7.34). 


Therefore, 


pe Do (a — k)* (4) 
k=0 4 
2 
ja7nial 
am Ly (m—k)?(,," 2) 
k=0 


Assume that the probability of a bit in error is p,. So, the average time period 
between two bits in error equals a In practice p. < 107%. As a result of this the 
expected average time between two bits in error is > 10007,. 

From previous discussions, we know that the effect of a bit that has passed 
through register R,; is the same as adding dc voltage to ie output. As a result of 
this and because the average time period between two bits in error is greater than 
10007, an ac coupling in the output of the RDC system will tend to cancel the 


effect in the output of the bits in error. 


E. SIMULATION RESULTS 
In this section we present simulation results when the RDC system is used to 
implement a chirp matched filter and a linear filter (lowpass filter) when additive 


noise is part of the input. 
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1. RDC Chirp Matched Filter 


The system simulated is shown in block diagram form in Fig. 7.12. The 
input noise to the lowpass filter is white noise with zero mean and standard de- 
viation sd,;. The bandwidth of the LPF equals that of the chirp signal, and the 


output band limited noise of the LPF has a standard deviation sd2, where 


3. z(t) = Asin (AE | 


> = o-0MzZ. 7 — 0 seconds.) — 8, and A = 1.0 





(7.36) 


When the band limited noise has a standard deviation sd2 equal to 0.75, the average 
results of six computer runs (where each run had a different seed for the random 
number generator) are as follows. The maximum value of correlation = 1.641; the 
output noise has approximately zero mean and a standard deviation 0.2768 where 
the average chirp signal power numerically is 0.472. See Figures 7.13 and 7.14 for 
the input and the result of one run. From the given results 


0.472 





(SNR): = apy = 0.762 dB 
2 
7 (SNR), = (a) = 15.459dB 


Therefore, the signal to noise ratio improvement SNRI = 15.459 — (—0.762) = 
16.221 dB. 

Fig. 7.15 shows the SN RI for different values of sdz; where theoretically 
the SNRI = Ty-B, and where Tj; - B is the pulse duration-bandwidth product 
of the chirp pulse. (See Horrigan [11].) Hence, for this example, the theoretical 
SNRI = 4(5) => 13dB. So, SN RI using the RDC method exceeds the maximum 
theoretical value by about 3 dB for this case. 

Denyer [37] introduced a digital technique called Deltic processing to im- 


plement chirp matched filters. We found that the SNRI using the Deltic processing 


i 





white noise 


x(t) 


chirp signal 


band limited noise 


Fig. 7.12. Block diagram of the RDC chirp matched filter that was simulated. 
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technique to implement a chirp MF is less than that obtained by the RDC method 
by 9 dB. 

The noise performance of the RDC system can be explained as follows. In 
fact, there exists a double match. The first is due to the matching of the filter to 
the signal. The second is because the delta modulator step size is adjusted to fit 
the input signal so there is no slope overload. When noise is added to the signal, 
sine sueiona occurs and the delta modulator does not always track the ine 


plus the noise. This lessens the effect of the noise input. 


2. Lowpass Filter 
In Fig. 7.16, assume that the LPF's have a cut-off frequency f, equal to 


5Hz and that f; equals 2.0 Hz where the input is 
i(t) = x(t) + n(t) 


where n(t) is zero mean white Gaussian noise with variance = (0.5)?. Assume also 
that only 5 zeroes at each side of the main lobe of the impulse response of the 


filters are considered. | 

Figures 7.17, 7.18, 7.19, and 7.20 show ive input z(t), and the outputs 
yi(t), yo(t), and y3(t). Fig. 7.18 shows that the DM acts like a filter because the 
step size for the delta modulation process is adjusted to fit the signal. The slope 
factor is s = 1.0. When noise plus signal are applied to the input of the delta 
modulator, the slope factor s decreases leading to slope overload distortion. This 
actually reduces the noise effect. In fact, the effective noise power into ie RDC 
system is less than the input power as derived in Equation (7.29). 


We can compare the noise performance of the RDC filter with a conven- 


tional LPF by defining 
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delta 
demodulator 4 
x(t) =1.0 sin(2T1 ft) yA ) 
- delta RDC 
= \ modulator ERs & y (t) 
2 
168) 
conventional 
3 


n(t) 


Fig. 7.16. Block diagram of the system used to simulate the outputs of 
the delta demodulator, the RDC LPF, and the 
conventional LPF. 
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dy (t) = yo(Ti, +t) — 2(Z) 
and 


d3(t) = ys(t) — x(t) 


where T;, = 1 second is the transient time needed to fill the n, stages of the register 
R,, of the RDC system. 

The simulation results show that the average power of d2(t) is equal to 
0.01324 where that of d3(t) is equal to 0.01902 which means that the output of the 
RDC LPF is closer to the original signal than the output of the conventional LPF. 
We then conclude that the “effective” noise into the RDC system is less than the 
input noise of the conventional LPF. 

Another way to show that the delta modulator is a “slope” filter is to use 


the block diagram given in Fig. 7.16 to find r1,r2, and rz versus o? where 


oF 
is anny 
O; 
a5 
rT = ~FZ 
oer 
o% 
rs = ay 
OC; 


and where o7,0%,0%, and a? are the variances of y1(t), yo(t), y3(t) and z(t) respec- 
tively. We consider the case where the input 7(t) is white Gaussian noise only and 
where the step size of the delta modulator is chosen to fit a sinusoid of amplitude 
1 volt and frequency equal to the bandwidth of the LPF which is assumed to be 5 
Hz. 

Simulation results are summarized in Fig. 7.21 as plots of r1,r2, and r3 


versus o?. Note 
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e For the RDC LPF, r; and r2 decrease as the variance of the input increases. 
This is because an increase in the variance of the input creates additional slope 
overload conditions. The delta modulator then acts as a “slope” filter. The 


slope of the output of the delta demodulator increases but at a smaller rate. 


2 


As a result of this the values of r; and rg decrease as oj increases. 


e When o? is small compared to the step size, granular noise occurs and so r; 
will be greater than r3. | ‘ 

The simulation results verify that the noise performance of the RDC sys- 

tem is better than that expected from the theory of linear filters. The reason for 


this is that the DM portion of the RDC system filters on the basis of slope as well 
as frequency of the applied voltage. 
F. USE OF VARIABLE STEP SIZE WITH THE RDC METHOD TO 


OVERCOME GRANULAR NOISE OF THE IMPULSE 
RESPONSE REPRESENTATION 


Aeeume that we want to represent the impulse response h(t) for a LPF or a 
band suppression filter using DM where the step size Ah is constant. For a given 
sampling rate, if Ah is appropriate for the main lobe (to prevent slope overload) 
then the DM values for all side lobes whose amplitudes are less than the step size 
are successive values of +1 and —1. This means that the side lobes are represented 
as being constant. On the other hand, if we choose the step size Ah appropriate 
for the side lobes, the DM cannot track the main lobe without increasing the 
sampling frequency. This See an increase in the number of stages of the two 
shift registers R, and R,. One solution is to use different step- sizes for different 


regions of h(t). We call this the variable step size RDC method. 
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1. Implementation of the Variable Step Size RDC Method 
To implement the variable step size RDC method, we first define A(t) in 


the range +47, seconds. To make the impulse response causal, shift the time axis 
by =+T;, seconds. An example is shown in Fig. 7.22. 

Now divide the causal impulse response into M different time intervals so 
that in each interval a single step size is appropriate. In Fig. 7.22 M = 5 and, 
because of symmetry, the step size for region A equals that of E and the step S1Ze 
for region B equals that of D. Assume that the M intervals have nj,n2,...,nu 
samples and step sizes Ah;, Ah2,..., Ahy, respectively. 


Using the procedure given in Appendices A and B, it is easy to show that 
BD Comings 15) = 7 eee ee (2 elie A hy 
me 
+ Sy [< Ra(t) >- < R,(2) >] - Aho 
meen as 


+ So [< R(t) >-< Re(i) >]-Ahs 


t=n1i+n2 
Bo eee 
SF eins. sus os + ven lguajecd « « 0700s 
N-1 
+ S > [< Ri(i) >- < R(t) >]- Aha 
w= N—n avy 
—eit R,(N) > is bas - Ts (7.37) 
where 
M 
N= a: 
and 
n;—1 nytno-l N-1 
is = ND) Ge aN: a es ra +Ahm: SA; 
i=0 i=n, i= N—-nvy 


= constant for a given impulse response and given step sizes. 
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Equation (7.37) can be implemented as shown in Fig. 7.23 where the constant 
Az -T, can be off the IC chip. 
2. Applications 

This section presents simulation results obtained using the variable step 
size RDC method to implement a lowpass filter and a band suppression (notch) 
filter. 

a; Variable Step Size RDC Lowpass Filter 7 , = 

Assume that the cutoff frequency of the LPF is 1.0 Hz. Figures 7.24, 

7.25 and 7.26 show the impulse response and the staircase representations (h,(t) 
and h,.:(t)) when we use variable and constant step sizes and where we take into 
consideration 40 zeroes on each side of the main lobe. The number of intervals is 
M = 21. Fig. 7.27 shows the resulting transfer function when we use variable step 
size. Also shown is the transfer function when a constant step size is used. Using 
the variable step size RDC method, the ripple in the transfer function is reduced 
and the filter skirts are steeper. 

b. Variable Step Size RDC Band Suppression (Notch) Filter 


Fig. 7.28 shows the transfer function of a band suppression filter where 


fp = 1.0Hz, jf, = Ss = 0.5, and Af = 0.375. 


and where 
H,(f) = Hilf) — Ho(f) 

Figures 7.29, 7.30 and 7.31 show the impulse response and the corresponding stair- 
case representation when we use the variable step size and constant step size RDC 
method. The number of intervals for the variable step size is MM = 17. The transfer 
functions obtained are shown in Fig. 7.32. 

From Figures 7.27 and 7.32, we can say that the variable step size 
RDC method has a better transfer function than that of the constant step size 
RDC method. 
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Fig. 7.23. Diagram of the variable step size RDC system. 


144 





e@eaenananee 


ececaaveaae 


dbeeeeeeeceeece 


aeseeteereeetdoserecececseeobvoseseasacashbsascavcvenectiaccecneenecececetbdacnan2eccaeccebaaeseeecneccetbaaaecetsearve 


Pre ccecce- 


eaetseeaneaene 


eeoeaaeareann 


weno regress consece 





OOS OF OS OF 6 0 OG HS OE we OOS OO FF OO OOOO OOOO FO SFOS ESTES |_ESSEF SESS SHES SHAMS STC SHEETS SSC CH SE BBEeTSSCeB etEEC2 


eaecaaeaeaee 


eeseaatoace 


eee ee2eqeeeeeeweeee 


earaant®@eaeene 


eaeeaenateeee 


eeeeeenaeeea 


ete eeeeses see Feoveeeceoceceteeeetosesesaetes este seeseeotesaestesetgesesseseSeeeeeeeceeneneo@®@oeseseeeceae se tevegaeeaasas 


o°O¢ 


‘gsuodsoal 


eaaaeseoerne 


*#eeeaneees 


eeecer daaanecasraceve 


eeeeea2rcaeace 


s72 a 


eeaeeeaeeaae 


eeoereezeeeeeeaaser 


eeeeece test agqgeses eer scarey esasetoceseqeseseataas 


asmnduit yesned Fd] eZ 


eaaaseaeene 


eeearnaeeoaoae 


ease eeeneececemes ee een eesn 


eeceecetoceyeretecatecsezreeeeeseceazosesescccecs 


eoearenasteen 


Coaaste F-feageaeh Fone art ems aeaeeeean 





er 


eee -eHRaeanGa ge Geo aosPaeetaanyaaeeate aes 


eee eer aa Gage fama aanrctamecaneanaee 


eeeteoacath yeh eoaneseetamereanetananee 


eaeeeseee 


efeeoereeeanea 


@eareeneee§ tC ge taoneeeceaemeeeneaeeense 


seers eee ee eteee eee este S*e at oatoer ars e®GaeeBarc@estera@esteastee eenteaneaeeneeaeneece 


eeeeenenenae 


SSCS SSS SS SS ESSCSS SS Sse Hse Fes seseseesessesespeteossesesesesgeoseseseseseeqeseseseneerqyseseeeneeeeaqgesece 


ecceccnccnc ze cccccceces 


ecccccewsaczocccccccave 


A 8 eg ee a 


N 
On 


qr eteeceresereqessesessovseyeseseeres 


; 


eresrecesoucy 


POL BY 


* 
* 
* 
e 

° 
e 
° 
e 
® 
. 


—_ 

ile me 

P ~~ 

Pee. 
CO 

fee By 

ae 

wn 
E\p oo cee 


eaerceateaenenaneaeeeeea 


= 
oO 
> 


. 
8 
* 
° 
° 
. 
° 
° 
. 
. 


° 
° 
8 
. 
. 
® 
° 
. 
. 
. 


® 
a 
e 
e 
e 
. 
° 
° 
. 
. 


e@eeeeeeene 


Sees eee ee eG ef 22 e2eeee2ees See eeeeet ett eet Pete aessteeset 2G atest eGeeeetateea* eb Bateaanearsese* ee eaveaanaacnea . 


+ 


ia 
* 
{ 


eaece 


@eeae 


eecacea 


e@eannee 


eteeenen 


eoeer esc ac etdaoaaneceeeetaaebeseeeeecevsesbasecacacesesastaaeaaservesbeeanesasraeatedeaacecesnasrecatesseecnesease 


NS 


ceeeeeenaataa 


@eoaoeeeeeeaasseeaeatonaae 


eaeeenneese 


eeeneaeeteon 


eaeeeveeeeaesne 


*' 


oT 


eaeeeceececeeeneeeeteaeee2eaneer 


eeeeeenececae 


' 
® 
° 
° 
° 
° 
* 
. 
. 
. 
. 


aeeeanereaanta 


eseereratgeeaeaatasene 





OOO 0 8 0 OOS OO RESO O SORE H|®SSHSSHHEHHSEHSSHSSSSHHSHHH|_ES SS SSSHT OH HEHEHC HT Hee FF eset EBearvrset** Fert gee eeeBBE7E8 


Paataancaeaa 


eaeaeteacacseosn 


Se 0S NS 1h 


ac UPC" U 


0S°0 


00° T 


S2*t 


+) 


14 








eeecececeeceetocaeeceve 


Ad] 24} Jo ssuodsar asjnduit yesnes ayy Jo uolyeywlasaidal osevoireys aT, “eZ ) “BI 


eaeeetotsteoce 





oe 2 





eect oeeeteageeree ee er eecepeseeeeeeene 





eeteeeececeoen 


eeeeeeceece 


eonrnervresveee 


eeoeseeereee 


ereetecrecee 


eotweraeestetacceeeesetcvoveveeetsnasasave 


eeeeeeteotavoeseneataoanetoantsa eran evant eonaaecas 


erteeeere008 


eeeeeeeecece 


TOTS eee essen yeresceonvys osc esseeyreeresece 


eerersreeeqrerecece 


ee) 
Ww 
€N}F oo ee 





eeceeteoceen 


eeerse tater 


eeeevetoecee 


e 
. 
e 
. 
e 
e 
e 
e 
o) 
e 
° 
1 
. 
a 
e 
° 
e 
) 
e 
e 
e 
e 


eee ee eee 


I 


eeeeneecerve 











eeoeeveeteone 


eoete®oaatoe 


eeeeeeeaees 


LIKE SO SMH I I MI HR HOLA St SKM A HEI ICN IK IM IID III IIE OC ICICI IOC IO OI OGIO OIG OO TON 
eee eeees OFee eee eee eFtoecesFte Sta FPeGeGoee Gee FFGusetoncvessesc e®ee esse seetoesesseeo999 F228 FFe8 298 


eneeeceeees 


eeeeteoeasease 


eeeeeenannee 


‘padojdura st azts days ajqeirea ay} UoyM 


eeeeaneton 


woeeeeccece Se eracve ees wecesee teat avec eeneetascecaece ete aneenves erst avesevee 


@eereeerececeegureveeceece 

) 
C) 
e 
e 
Co) 
° 
s 
4 
e 
4 
« 
. 
s 
. 
° 
) 
) 
s 
e 
e 
° 
a 
s 
. 
e 
s 
e 
s 
® 
* 
° 
° 
° 
1 
s 
. 
° 
e 
s 
e 
s 
e 
. 
. 


ereeeseeeeesmeseseseanete 


: 


eee4ee®eeneee0 


eeeeteseeeen 





eereteneeee 


eeeteoaaeone 


@eteoenaete 


erceoeeeeeoeemgeeeeeeeeteosr 


eeeceseeeene 


esreeeaeae 











eeerseert pean pee Meeeeeass 


eeeeqe ene tee tegen eseesece 


eocee G-eoeohoto hae MH eoeee 


efe we eeaeeete 


eaeeeetresk pgeaneeaeees 


eeeeeeeesetaoeeeee ee foes ese veeteseevenes Peano 2988 





@P@aereaeeen 


OS OS OS Oe 98 8 OF FOO SETS SSS THT ET SHHSSOTS SEF SOSH SOS+SE STFS SSFSSSS|EGCSSSTSFCSFS*SSFTTCTCFCFTSRETVISSS 


eeaeetoeoonen 


eeoaeneeeonns 


eerecsceenece 


eeeeeeeeee 


eeereerear sth eeeeteee ee mete eceseanenen 


eteeaen te 


eeeee0e0e0 


Tt HR th tS I I TIS, 


‘ 


Seer ev ese_zveserr eae 


eeeeccseegy 


@eveeeenne 


@eaeeeenoes 


S@Geeeeeecece 


at core Se at oli 


es 
e 













0O°t S20 08 @ Sc 0) U0 ORs¢ 005 te - 


fe t] 

i) * 

8 ¢ 

° e 

° ° 

° ° 

° 8 

* + t] 

8 i 

eeeC¢pPeeeteae e 

° CY) 

° e 

° 8 

e 8 

e t) 

® t] 

° i 

° 8 

i) t] 

BPeeeettoeceeaeeeteereaee S 

e 8 

e ° 

® e 

° 8 

8 8 

s * 

° t] 

t 8 

e + 

er2eeeeacece ° 

° 8 4 8 

e 8 8 8 

e J 8 e 

® 8 ° ¢ 

4 CY) ® ‘ 

® e ® ° 

° 8 ® iy 

° 8 ® 8 

° 8 . 8 

Geeeeeneteageranereecezecesn Raeeonseereerace i @reeetoveets 

° 8 8 i] e 

° . . 8 ° 

i} 8 8 ¢ ° 

° 8 e ty ° 

° 8 8 ® e 

° 8 8 ° 6 

° 8 ° ® iy 

e ® e ‘ ‘ -—> 

° ° 8 . e 

@eC¢evpeeeaeeven @eeeeereeene meeeeGeeace @ ere aoaeeaeeete @r®seeeeteee e 

° 8 i 8 6 

: : : ° : 0) 

e 8 % + iJ 

: : : : ‘ on 

e 8 8 ' CY) 

° 8 8 ® i] 

° 8 8 8 8 

® ® t) 8 ® HD 

e 8 8 ° s 

eePpeeeeaeresn Peeeneeeere Reaeteeseanne wee Peeeteonaee Ser®eee@eeee0 e 

e 8 8 8 e 

° ’ ’ : : SA) 

° 8 t] i] s 

° : : : : O 

° t] ° e e 

® ¢ 8 ® ry 

e 8 8 ® ° 

® 8 e 8 LY) D> 

° ° 8 8 i 

ere eeeneete Seeececceccce Reet Goeeeee ePeeteaesece @weceeerenee e 

° 8 8 ® 6 

® 8 8 e e “SJ 

® e ® ¢ t] 

° 8 t} ® ty uw 

° i] 8 ty 8 

® 8 8 ' s 

° i] t] ° 8 

° 8 ® ® 8 

: : , ° : N 

@eSCSaGaeeeea eeeeeeste eee gteteean Set eee oi eee Saoaeneteoat @eeeeeeseten e 

° * 8 8 e 

‘ ° ° : ‘ © 

® t) ] e 8 Ee 

® i] s i] 8 

L} 8 8 ® 8 ray 

e . ty ° 8 

° t) C) ) 8 ( + ) U 

® i] i] 8 8 

® i] 8 8 ® N) 
es 


Ge 


6 


14 





eee er%® eeoeae 


Cee ree eerere ree eres) 


weaeeeessesatdos cart socateboseeeeeeesebocseeeeesectawassemensetdaeeeeseseeebdesreeaaeeseesasesese 


‘poAojdura st azts days yueysuo. ayy Us M 


Ad 243 jo gsuodsal asjndult yesned at} jo uoryeyuasaidal aseoire4ys 


eeeer res ee op seers eneene 


wee ereeeeeeeqgeeerees 


eer eet eet otpeseeevneseeen 


eeeere2aven 


e@eeeeeerece 


ee eeetoeeqgeseseneaeerece 


eereeceececeregaeenreceeaeyene 


eeeneeeens 


SOS SSC Sse FO Yss Fees SH sOS Goat atest ests ao geet sts ost stat get eaese* eo eeegseseeeeeeeees_eeeeseeen 2eeqeaznare 


8 

? 

8 

ry 

. . e ‘ é 

e ‘ ’ r) \B ° ‘ 
ry CY) ‘ , . ‘ ° ) 
e ‘ » ¢ r) ‘ e 8 
C) ¢ ‘ o] e ‘ e 6 
e C) Cy) e t) 8 ° e 
° CY) ‘ ‘ t) CY) ° e 
eeceervreeee VPP Pm et esta Best Foe eeeeeee ees Sa Faapeca- eee Cot ammteose tse Faetpeneeet ent ag 
6 . ‘ ° . ¢ ° ‘ 
e ) ) ? 8 . ° ° 
8 ) 8 e t) t) ° ° 
‘ e ) t) ? e ° e 
e ® CY) e ’ ° ? ? 
C) t) a ? 8 t] ? ? 
‘ 8 ‘ r) CY) 8 ° ‘ 
e ® a r) 8 8 e ° 
e 8 t) e r) 8 ° ? 
e e ‘ ) 8 ® ° 8 
e CY) ® e ty) C) ° ) 
@eseeeeeneee e@ereeeteeenee Q®®eeeeeeeerpeseeeeteavergeereateaane ee? e@eeoaneeaeemszr eee eee ee etempaeteaaeeeeteg 
e i} ry ° a * ? i] 
° ) 8 ° t) CY) ° e 
e . ’ e C) C) ° ° 
e ) ‘ ° r) CY) ° ‘ 
e t e ? 8 ? ? e 
8 8 8 ° ° s ° ? 
ry t J 8 L} e ? e 
° ° ty) ’ ® r) ° 8 
e ) CY) e r) 8 ° e 
e CY) 8 ‘ ty) t) ° 6 
r 8 ‘ C) t) CY) ’ e 
eeeees Stes egrets sGonee QF Sesto auetpoeseraseenepxPeeeeeanee Te gtx a serateceqgrtosnetoeneeecempeaseeee2es vg 
C) CY] 8 ' CY) * e C 
CY) t) ¢ ¢ 8 8 ? ? 
r) C) 8 ? 8 CY) e ° 
8 CY) ty) ° ° ’ e ° 
e t) 8 ’ 8 t) e ‘ 
C) e ® ' 8 8 e ¢ 
6 CY) 8 t) 8 t) ° 8 
e t) ‘ s ‘ ‘ ° ° 
e CY) CY) ° CY) CY) ° ’ 
) CY) 8 ’ ‘ C) ° ? 
’ ) CY) ’ * ’ e e 
Pete etee see Geee ese Geos eae eeeeeeanat ewes ee Gene teGSeet®eannt®s *h eo Geaneennevace ~s~eeepeeeteeaaetFeeteoevnen et oop 
° . ty) ' ? t) ? ° 
° CY) CY) e ° ? ° ° 
° 8 ® e 8 . ? ' 
° ® ® ° e e ? ? 
? ? ° s 8 ? ° ° 
’ Ct) e ° 8 ? ° ° 
r ) ® ? ? ) e ° 
? ) LY ’ ‘ ° ® ° 
t t) ) ° ° * e ° 
* 8 t) ° ° 8 e ° 
? 8 8 ’ 4 e ? e 
e@eeeeevaneaene @eeeereeee eae Ber eee anne aeer ete ©* Gea eeaarGeeoeenvreaanenanse e eeeeteeeoeePmeeaeveeeneeeeenemeseeaneeeseee 
? CY) ® ° ° CY) e ° 
’ t) 8 ' . ° Y) ’ 
e CY) ® e ‘ e ° 8 
¢ t) LY ’ 6 ° ? t] 
? . 8 ? ’ e ? ¢ 
‘ Cy) r) ? ‘ ° e ° 
’ t) ® ’ ‘ ° ? ? 
° C) CY) ° ° 7 ° ° 
° ) 8 e ‘ ° ? ° 
? ) ) e e ? ct] ) 
6 CY) 8 e ° ° 8 8 
eanneeve S@enteseeerat eat eReeneeeeeeears ee eeeeeeeersePeeeenenatves @eeeneeeeeme nen et et eveeFeae sees eearve op 
e t) ¢ r) ? ° ) 
¢ t] a ' ) ’ ) 
¢ ty) ) C) ° e CY) 
¢ ty ‘ i 8 ¢ e 
° ) e 6 e ° s 
? . e 6 ? ? 8 
> e cy 8 * 2 ? 
* ) ' 6 8 ? ? 
° ) . 8 ? ° ° 
e ) 8 s ' e ? 
e e ) ' ? ° . 
eeeqeeoenestege ese@®@eeeene Qe re @eeeavepeteteoeveeaneanse pPeeettoee Gata k wee aeneees& meeereeev ene opmeecenese @erveeg 
° ) . CY ) e ' ° 
i! e ® C) CY) ? e . 
r] t) . ) ° ? ° 
? CY) ’ . ‘ e ° 
° CY) ? 6 e ? ’ 
° ) ° 8 e ) . 
° ® e 8 . ° . 
? t ? ? e ? ° 
° s ? ’ ) e * 
? ® ? ) ? ¢ ? 
? CY] ? ? ? * ® 








eeeeseaeeee Q®*FSeeeaeeeteopreereeeeveegeeervrae § Bea ga Of Foe eanes 





ee eee eee ee 
eo eo es 
ere ewe greeter eaeeae 

ec swedtdwonon sc eeee 
Ce Ce 
eeseeeegeceeeecenesen 





(2 





mae. 
=) 






eeeetetvertece 





UL WL 3 


eeeceeee ee eee 2 eeee2 eee 220848 SB we em ase eeee2e0 


ee®esreeeeneebececenesnceeastoseeeeeesesets eo cases 


eeceeceotasaneereseoee 





Mieeerser22 een et eeee8 


eeoe 


eee @ereeeten 


Seer ee ee Seaegr ers ae Sr ses ert ee eeerst ester sagteeuce 


° 
) 
® 
. 
. 
° 
e 
. 
o 
° 
ty 
e 
. 
ty 
a 
) 
® 
. 
e 
e 
* 


eeee 


®@@eereeneeeceeeeee eee ee wees eerteaeeearsee eee eeeen eee eeeeeseBosn 


LI 
e 
. 
° 
. 
. 
. 
° 
. 
e 


earteeeteoce 


Terre rereerrrrere ere rere ee re eee eee eee ree eee ree ee 


eese seen e®? a 


grec eeesrace 


Seeterene=ee 


seeevreevreans 


eee ertaerteee ger 2 Fs est eFest Freee Seseeenegeeeeesese een 


Se O- 0S 0- 


Geel) y O07 ol Cee use0 Se 0 00°50 


OS*] 


00°24 GZ°T 


Ve 


147 


ececereeeeeeeen 


¢ 


” 
pote w ce ws ecw cee eee-- 


e®eeet®t®eeaveat 


eeeeeesreereeece 


Setewmese ese ee eet e Peet eee seee sees e eee C2 FF ssF ees FFF est soMeoesetsessFeoeotstessessteoFseeteseezesenn 


(at a: 
Qs 


eesevreretee oa 


0] 
v= 
N 


= 
ee 
e 


‘STI OCU 221s days ajqetrea pure yuRysuOD jo suolyoUNy Jojsuety, “LoL ‘Sig 


. 
eoo 74 eeeerseee 


e 
e 
e 
es 
e 
s 
e 
8 
° 
+ 
s 
8 
e 


e 
. 
¢ 
¢ 
e 
e 
° 
e 
e 
* 
. 
8 
. 


eee eeeeteeaeat eevee avraeneaartresctseeseaaeaeevreartePaeseGeaseeseeeeGatst at eavetstetvaasi®tersereeras eee ee PGeaveseeneneece 


Fa eS IIS 9 EN ST SC) 


° 
6 
s 
¢ 
e 
e 
se 
s 
s 
* 
s 
i} 
* 
. 


eoeersroeotoeente te tease eetevreeeFevteeneaaeavrevrbeseeaeneaanavrees 


a 
° 
) 
. 
s 
+ 
° 
e 
° 
» 
. 
® 
° 


eoeeeeeeeeeceePeaaet eee eeave 


° 
Terr rrerrrre eae Serre eee eee eT eee) 


oo” 


CIS IT SII BI II OOK O OS HII ITI OCICS I CY To ST 
-— 
e 


eeeereeseeecsceeehb eee e pee et eee OF eOeeeeseeeeeseeeetaereseeeet ee Fase aseereeereas Fen seaseaeeeeee Ge aetoevaveese seh eae deat eeeeeon 


8 
e 
° 
e 
. 
s 
8 
¢ 
s 
s 
s 
s 
. 


seeeeee eee Peano teevee®obBeeeeesn 


eee ee eee ee eet eee are FSG et eeeeeeeeeeteneeeneeeeeeeeeasen? 


eeeeeeeee emer r er ee eeeaee ee aeeeaeaneereenaste 


wo ; 


” 


«ese 
Cie 
- 
° 


1 
) 
s 
s 
. 
, 
* 
Co) 
* 
+ 
CY) 
s 
, 
eeeeeee Meee eeeeeeeeeesceeweee st eee est F FFs est FeMesee Stet eereeetenehoeseseveecenzn ee eesesbesseereeseeeneees 
. 
° 
s 
s 
) 
s 
s 
s 
s 
s 
» 
6 
, 
Sere Fe oe PP ss ser Fse2FSessOoM er eee sv es ess eves eH eetestevevesseseseseFseerre Feet ee Fees stsseeresesevetstesesss 
° 
s 
s 
e 
s 
> 
s 
* 
e 
) 
s 
° 
. 


* AN 


weeeeeeececes 
° 
- 


ceeoeoee geese e@e¢ee @eQeest seer eevea eeeeeeeteeece eeeeeeeeeeaete 


SE gALS INBISNOS 


eeedseaee 
e 
e 
» 
°* 
s 
e 
° 
s 
2 
° 
° 
“ 
¢ 
° 
e 
s 
e 
e 
s 
°e 
s 
s 
e 
e 





weteteereos ee eaeee 


eevee ee ees etM#eseee 
« 


Siete Ltaataieiare 


s 
e 
® 
‘ 
s 
1 


ee ee ee eee 
—_—oO = 


gc mg 


"7 


em eee ertoeee 


weowe eee ee eeese 
® 
s 
° 
. 
oi 
® 
* 
° 
° 
, 
el 
. 
ef 
. 
l 


=) 
o 
(=) 
oS 


weeeeeeceeee eer eewebe see eeseeeeteneviecaneess essere sees 


. 


~ 


®ereeeeeeteaneé 


se eereeeeeeeeee seer ee eee eeeene eee eee tear teygteeeereeeetaevegzetoeaeetheaseeeetpPeseeFesrees ese eepeseeeeeeeeeeer~s Fever eeeeesee geese @eveeeeere ye eeeeseeneouave 
A 
n' 
s 


ere teteeeeen 





0° 52 


0 Sc- 0°CL- 


0°0¢- 


0°S- 


(dP ) 
(F)H 


148 





Fig. 7.28 A band suppression filter transfer function H,(f) where 
Hy(f) = Hi(f) — Ao(f) 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The RDC system is a new realization of digital filters. The system can be 
used as a convolver or as a correlator. The RDC system has a simple layout for IC 
ealinateon When compared with other digital filter realizations, the RDC system 
appears to require less hardware and fewer connections. The system requires no 
synchronization, has no multipliers and is easily programmable. The RDC system 
provides real time operation (no off-line processing). There is no limitation on the 
time duration of the signal to be represented in matched filter applications. Simu- 
lation results indicate the RDC system achieves noise performance which surpasses 


that predicted by the theory of linear filters. 


B. RECOMMENDATIONS 

It is recommended that a RDC system be fabricated and tested to confirm 
simulation results and to provide additional insight into this method of achieving 
convolution. 

It is also recommended that additional study be done and the results be ex- 
tended to include adaptive delta modulation techniques. | 

Also, studies of other noise reduction methods including additional work on 


slope filtering are recommended. 
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APPENDIX A 


IMPLEMENTATION OF A TRANSVERSAL FILTER 
USING DELTA CONVOLUTION 


Assume that we want to convolve the analog pulse signals A(t) and z(t) where 
h(t) and et) are the delta modulation versions of A(t) and x(t). Assume noes h(t) 


and z(t) eve number of bits n; and nz respectively where 


ny <N2 


Assume also that we have two registers R, and R, where FR, holds the delta 
modulation bits of A(t), and R, holds the delta modulation bits of z(t) which have 
to be convolved with the contents of R,-So, Rp, has n; memory stages and R, has 
Nn stages. 


By expanding Equation (3.18) we can write the result of convolution at time 








7° as 
Se! J 7= hy 
PesaG. Le) =Car h at dat das + Ya 
0 0 7=0 
j-2 j-3 pny 
+hy Soait Soa + tq py =: 
2=0 1=0 2=0 
g=3 aed 
oe Sees c 
7=0 2=0 
sha! mails (acest M+ gens age 
Se a es 
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which can be rewritten as 


y(t) ~~) yQ ° Te} Oe fc ej ap + JT ;=3 ee ee 


+hy | Lj—2 + Cis Se eo 0. bus,’ 
= j-™o 
+ (ny —1)- 27m, + (m1 = 1) x 
=) 
SI ae fre eh, Co oe 
5 ee re en cue. 5 oe. 
j—ni-l 
+hn, —1 by =| 
1=0 
(A.1) 


The products needed to implement Equation (A.1) can be represented as 
shown in Fig. A.1. A schematic diagram of a system which provides y(j - T;) 
according to Equation (A.1) is shown in Fig. A.2 where f, is the sampling fre- 
quency. 

In Fig. A.2 we assume that AR, register stages can have any of three states 
—1,0, or +1 and that it is bidirectional. Assume also that we have a logic which 
can do the following: | 

a) During the time periods between even and odd sampling pulses the contents 
of register R;, move to the right. 

b) During the time periods between odd and even sampling pulses the contents 
of register Rp, move to the left. 

Register R, has 2n, stages. The first n, stages have to be initialized with the 
delta modulation version of h(t) which can be either +1 or —1. The last n; stages 


of register Rp, have to be initialized with zeroes. 
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pf 
DD 
y. 


Fig. A.l. Diagram showing the products needed for delta convolution. 
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From previous assumptions, during each time period between even and odd 
sampling pulses the contents of Ry, move to the right n; times producing all prod- 
ucts needed to find the convolution, and during each time period between odd 
and even sampling pulses the contents of R,_, move to the left n, times producing 
all products needed to find the convolution. The contents of the accumulator are 
multiplied by the constant C, to produce the convolution. 

Some of the disadvantages of this direct implementation of the delta convolu- 
tion are 

1) By looking at the columns of the XNOR gates of Fig. A.2, we see that 


the needed number n3 of binary multipliers (XNOR) gates) is 


nz = ng +(ng—1)+ (ne —2) +... + (ng — 7; +1) 





n2 n2—Nny, 
ci a: 
7— {= 
= (mi) 
2 2 


2) The memory required is 2n; + n2 bits. 
3)-The high clock rate n; - f, needed to clock Rp. 
It is clear that the required number of XNOR gates and the required memory 


will be large, especially when nz is large (which is a practical assumption). 
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APPENDIX B 
MODIFIED AND REDUCED DELTA CONVOLUTION METHODS 


1. Modified Delta Convolution 


We will derive the change in the discrete convolution Aconv(z) when the input 


_—~ 


register R, is fully occupied where 
Aconv(i) & conv(i- T,) — conv ( —1|- nd (Bas 
From Equation (3.14), at time t 


Wt) (kT) = (Gi +m)-7s) =C, : 


ho - ae + 22 j+n,—2 + 3Z j-+n,—3 See + (ny _- 2) “2542 (ny — i) 
j 
“Zp41 TN So zi] 
7=0 
Eon © | Pipa, —2 + 2554 my 8) tues Petit ea oe) 


j 
*2j41 + (my -—1)) 0a) 
2=0 
+ h2- | Ljtn,—-3 +--. + (ny — 4): tj42 + (nr — 3) 


A, 
-Zj41 + (ni — 2) 3 ra 


2=0 
Se a Ais os oe 6 Sc 
SREP ER ERO rs COE SAGO Ss ook oy 5 6c 
j 
ate eat >] 
7=0 
where 


: t 
k = integer value of =| =jtny 


3 
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and where the needed products can be demonstrated as given in Fig. B.1. Collect- 


ing terms gives 


ue) =U +m)-T) = 04-4 
pany —1 j-l 
(| 2 zi-(ny —7 +7) rhoonSa 
. j+n,-2 
(ns | 3 zi-(ny —7+97-1) 


7=0 


+ hig(my —2)- r | 


i= 





+hy- (ny —1)- =") 
jtn,—3 
+(m | S> ai: (my i +75 -2) 





=) 
Jue, ee 
5 = SR oe eR es ee 
j jul 
tha,-1- |S ai-(n1 —i +9 —(m1—-1)) thai Dai } 
i=j i=0 


n,—1 j+n,—m-1 
yk: T) = Cr » (1 | =. nto ij —m 


m=O =) 


- (S's) : "S. helm =) (B.2) 
7=0 L=0 


At time (k + 1)-T,, 
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Fig. B.1. Diagram showing the needed products to find the convolution 
at time (j7 +71) Ts. 
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y ({k-+ 1] Ts) = y (lr +i+UT)=Cr-4 








ho : Ea + TED Bc —1 + 3 j4n,-2 = oer + (n, = 2) "L543 
eed: 
+(m; —1)-2j42 +71 - ie 
i=0 
met | Ljtny—1 + 2Fj4nj-2 +... + (m1 — 3) - 2; 43 . 
pat 
+(n, —2)- 2342 +(m1 - 1): doz 
1=0 
+ h2-| Ljtn,-2+--.-+(n1 — 4) -2j43 
Sere 
or tee) bye + (7 — 2) - ee, 
sto Sel oo SUD One Seo 0 Oe Re 5 6 Oe nr 
jt+1 
ate ha, —1 | ya 
:=0 
(B.3) 
Therefore, 


y([k + 1] “73 = y([n +j+1]-T) =, : 


Jtn, j-1 


pent I=} 
ar (, : | sy ri-(ny -i+s) +h, (si +m -1-E2)] 





jtn,—2 j-1 

- m| » no -ij-0) +h (-2 +h 2 Ea) 
i—7 1=0 

ae OR ee oe ee 
gt+1 q-1 

+ | hn, -1 So aj- (ni i+ 7 — [ni —2]) + hy, -1 “+305 \ 
i=) 1=0 

Or, 
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ni—1 jtni—m 
(eta) =6-4 Ss (| ye nolo nitj—m-+) 


m=0 i=—o 


in (six) -(¥ busty =H) =) > mt (B.4) 
L=0 1=0 


1=0 


From Equations (B.2) and (B.4), now form Aconv(j + n; + 1) where: 


Aconv(j +n; +1) = conv({j +21 +1]-Ts) — conv ({j +7;]- Ts) 


n,i—1 jtni—m ni—1 
=c.-({ S (1 | > nolo iim) ee » st 


m=0 =] 
Nn; —1 jtny—m-—l1 
=) (he [SE situs its —-m)} P| 
m=0 i=j | 
ni—l jJtnim—m 
=c.-( te SS rie(ny —t+j—m+1) 
; m=O =) 
jtni—m-—l : ni—l . 
— > se-(m -itj—m)| baa a) 
— 7=0 
nil jtn,—m—l 
=61-(| 5° hn ( S° slo bj m1) (ma 845 my] 
m0 i= 


n1,—1 jtny, re ni—l 
0. (|5 m ( s nttiemen]| oa ay 
1=0 
: n1—l jtny—m nil 
«(ECE Es 
m=0 i=] 7=0 
n1—1 jtn,y—m 
-.(§ ma ( yy “)), (B.5) 
m=0 i=j+1 
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Or, 


Aconv(jy+ny +1) = Cy - ( yy <itp(m) >: ( ~~ a i) >) | (B.6) 
m=0 = 


where < R;(z) > means the contents of cell 2 of register / where k can be fA or vw. 

We can see that the change in the discrete convolution is dependent only on 
the last n; samples of the input and is independent of all other samples. Also we 
can see that all the weights equal unity. The needed binary products to find the 


change in the discrete convolution are represented in Fig. B.2. 


2. Implementation of a Transversal Filter Using the Modified Delta 


Convolution Method 
To implement a transversal filter using Modified Delta Convolution, two 
methods can be used. 
a. The First Method 
In the first method (Fig. B.3), the needed binary multipliers (XNOR 
gates) = ny + (nm; —1)+(n1 —2)+...... = = teeth = oieyots and the 
needed binary storage =2n. 
b. The Second Method 
In the second method (Fig. B.4) we will use the fact that we need to 
perform n,; binary multiplications onlv and to store all the needed binary multi- 
plications which were performed in the preceding steps. In this case the needed 
number of binary multipliers = n;, and the needed storage is as follows: 
2ny +(n; —1)+(n;, -—2)+...... +1l=n,+ ea — eu Despite the fact 
that in the second method the needed hardware is reduced, we still have to use a 
number of XNOR gates or a number of binary storage locations proportional to 


(n;)?. Because of this we introduce the Reduced Delta Convolution method. 
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rocccc- 


<R (2)>=x 
x 







Jeny-2 


oo lh es % sen = | 


Sle (n -4) >= X jet 


Fig. B.2. Diagram showing the needed products to find the change 
in the discrete convolution. 
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. 








Fig. B.3. Diagram showing the first method of implementing 
modified delta convolution. 
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Fig. B.4. Diagram showing the second method of implementing 
modified delta convolution. 
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3. The Reduced Delta Convolution (RDC) system 
In this section we introduce the reduced delta convolution (RDC) method 
where the needed number of binary multipliers equals the filter length n, and the 
needed binary storage equals 2n, + 1. 


Let us define A2conv(k - T;) as follows 
~- A2conv(k - T,) & Acon(k-T,) — Aconu([k—1]-T,) - {B.7) 


where Aconv(k - T;) is the change in the discrete convolution which is defined in 


Equation (B.1). 


We now find A2conv(k - T,) when the input register is fully occupied. At 


time t where 


t 
k = integer value of =| =jtny 


we know that 


Aconv(k-T;,) = Aconv({j +71] -Ts) = Ci - 


-— ho [zjen,—1 + Bjtn,—2 + Bjen,—-3 Heese: + 2;-1 +2;] 
sore Mee aor tee aut... + 23-1 +2;| 
Sey iy Ds. re a eae ia Z| 
=P obo oO Ne ba ond Sienna 
Foc 0 0 0 St ERE OG. 5 Ge 
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and 


Aconv ((k + 1]-T,) = Aconv({j +71 +1])-T,) =C, - 


ho [Zs4n, FSpeageiee oe eee ee +2; 4+ 2541] 
+h, | Ljtny—i or 24a 2 eee +2;4+ 2;41] 
+hz [ a Ee 0c +2; + 2541] : 
errr rrr oo no SU oo 
on cee 2s iiss i oe 
+hn,-1 [ snl | (B.9) 


From Equations (B.7), (B.8), and (B.9) we have 


A2conv ({k + 1]- T;) Neonat +1}]-T,) — Aconv(k- T;) 
A2conv ({j +n, +1)-T,) = Aconv ([j +n, +1])-T,) — Aconv ([j7 + 24] - T;) 
=C\- {ozsen Se ee rae OE aa oP ae oo 
ni —1 
1h aoe Ss | (B.10) 
7=0 
which can be represented by Fig. B.5. Referring to Fig. B.5 and Equation (B.10), 


then 


ni—l 


A2conv ([j +n, + 1]-T,) =C; - > ee eo Decay 3 
1=0 
ny —1 
a R,(n,) oe _ st (Bales 
2=0 
From Equation (B.1), 


Aconv ({j +n, +1]-T,) = conv ([j +1] -T,) + Aconv([j +n, +1)-T,) (8.12) 
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and from Equation (B.7), 


Aconv ([7 +n; +1]-T;,) = Aconv({j +1] - Ts) + A2conv ({[y + ry + 1)- Ts) 
(B.13) 
Se eee celts en by, Bonation (B.11) and where at 
fime? — 0, conv(0) = 0. 
Soe eds eae Mer AES enna wnat welcall the reduced delta 


convolution (RDC) system. See Fig. 4.2. 


jul 


=} 


<R (na -1) >=h. 





<R. (0) > 
XNOR gates 
from DOM 
Sees Xx ee, 


<R (0) >=x 7” 
~ : <A (n-1) >= ng 
See os _ 
2) 2h \ J 
+ - 
A, K A 
FromXNOR = + . : 
gates a 


A 2 conv(j+n pd) 


Fig. B.5. Diagram showing the creation of A2conv ([j +n; + 1] -T;). 
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APPENDIX C 
A CONVENTIONAL MATCHED FILTER 


Fig. C.1 shows a block diagram of a matched filter where 


Ta Z 
y(t) = / h(r)- w(t — r)dr 


0 


Ta 
ie z(Ty —7T)- w(t — 7)dr (Cal) 


0 
where z(t) is the signal to which the matched filter is matched and Ty is its time 
duration. 
The first approximation involves writing a sampled-data form of Equation 


(C.1). See Equation (3.15). 


L-1 
yy, ~ te SS TL—; ° Wk-i (Gz2)) 
1=0 


where 


Y; = yy ; Hs) Ob = w(J ; ae), and Lj = AG : ihe) 


The sampling rate = = must exceed the Nyquist rate. L is the number of samples 
~ and so Ty in Equation (C.1) equals L- T,. 
The samples z; and w; are now represented by binary words having elements 


1 or 0 identified by c® and w” where 


oo Stem 27m, w; = Swe 2" (C.3) 
m=0 n=0 


where 2 and w? € {0,1}. 
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w(t ) a) 





Fig. C.1. Block diagram of a matched filter. 
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We now approximate the sum in Equation (C.2) by truncating r; to M digits 


and w; to N digits. As a result of this 


M-—1 [N-1 L—1 
yk © Ts - Ss 3 eer) 3 (aig no olga (C4) 
m=0 Ln=0 2=0 


Correlating the m digits of the binary words of the samples of z(-) with 
the n#@ digits of the binary words of the samples of w(-) gives. the digital corre- 
lator (m,n) given in Fig. C.2. At each sample time the most recent sample’s 
nt digit enters the left most stage of the n on-line shift register. The m™ dig- 
its of the sample z; are permanently stored in the stages of ae reference shift 
register. But, each sample w; has M bits and each sample zr; has N bits. 5o, 
to obtain yz we have to perform M - N digital correlations, where the results 
of correlation are weighted by the factor T, -2~("+"). We use one reference 
shift register, which has M- DL binary stages, for all the M-N digital correla- 
tors. | 

Fig. C.3 shows a block diagram of the DMF for arbitrary N and L but with 
M equals to 1 because a block diagram for arbitrary M,N and L showing all inter- 
connections is cumbersome. In Fig. C.3, for convenience, the gain factor T; called 
for in Equation (5.8) has been left out of the weights 2—~"*"), so the output in 


ee Ue. 
fact is T° 
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eeees 


the nth on_line shift register 





reference shift register 


Fig. C.2. Digital correlator (m,n). 
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digital correlator (0 , 0) 
digtal correlator (0, 1 ) 








w(t) 


digital correlator (0 , n ) 
digital correlator (0 , N-1 ) 


Fig. C.3. Block diagram of a DMF for arbitrary N and L, where M = 1. 
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APPENDIX D 
DC OFFSET LEVEL 


After the input expires, the zero voltage input to the delta modulator will be 
represented by a successive positive and negative pulses about zero. From Equation 


(3.18) the output of the convolver equals 


n,—l jo—-k-1 
yi) ey Ti) =e aie 3 {| 3s) b> 2 j= 1, 2.35 


k=0 1=0 


But the input signal delta modulation pulses z; have successive positive and neg- 


ative values. Therefore, 


j-k—1 , 
3 Pe! ce for j —k = odd (D.1) 


for j—k= even 
7=0 


From Equation 3.11 


k 
h(r) wh(r)=Ah-S hp; &-T, <7 <(k+1)-T, (D.2) 


p=0 


Therefore, the output level at a given sample time 7 -T, when 7 is odd (k even) 


becomes 


Te ys h(k-T,) | -(+Az) 


k even 


At T, seconds later, 7 will be even (k odd) and so the dc offset value is 


nyr—l 


HO || SS Ae Ty || (eeu) 


k odd 
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Because T, is small when we use DM, then we can expect that the distortion 


will be constant and approximately equal to 


aa (total area under h(t)) CD73) 


since T,, - 3 h(k -T,) is approximately half the area of h(t). 
k odd ae even 


Lig 


lO 


MO) 


1 


Wey. 
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